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I.  SUMMARY 

The  seismic  motions  from  nonlinear  two-dimensional  (2D)  axisymmetric  finite 
difference  calculations  of  underground  explosions  are  propagated  to  regional  and 
teleseismic  distances.  Synthetic  seismograms  are  computed  and  examined  to  study  the 
effects  of  depth  of  burial,  medium,  structure,  and  prestress  on  seismic  waveforms. 
The  elastodynamic  representation  theorem  is  used  to  couple  the  results  of  finite 
difference  source  calculations  to  regional  Green’s  functions  computed  using  a 
wavenumber  integration  algorithm.  The  2D  finite  difference  source  calculations 
include  the  nonlinear  effects  of  gravity,  spall,  bulking,  and  cracking  as  well  as  the 
finite  extent  and  divergence  from  spherical  symmetry.  The  results  are  broadband  full 
wave  regional  synthetics  for  explosion  sources  with  cylindrical  symmetry.  These 
theoretical  seismograms  are  studied  with  regard  to  proposed  regional  discriminants 
between  explosions  and  earthquakes. 

Synthetics  for  Pahute  .Mesa  tuff  are  examined  for  three  depths;  nomtal  burial, 
overburied,  and  cratering.  Regional  synthetics  from  axisymmetric  explosions  are 
significantly  different  in  the  short  period  baud  (0.25  to  2  Hz)  from  simple  point  source 
explosion  models.  Simulation  of  an  overburied  explosion  shows  enhancement  of  high 
frequencies  relative  to  a  normal  buried  source.  The  2D  axisymmetric  sources  exhibit 
greater  Lg  excitation  than  point  sources,  and  the  normal  buried  source  shows 
enhanced  short  period  fundamental  excitation  compared  to  a  point  source  or  an  over¬ 
buried  source. 

Regional  synthetics  for  a  simple  model  of  Shagan  River  explosions  also  demon¬ 
strate  that  the  point  source  is  an  inadequate  model  for  regional  phase  excitation.  By 
contrast,  regional  synthetics  for  model  explosions  in  granite  with  or  without  prestress 
differ  little  from  the  simple  point  source  model  in  the  short  period  band.  The  implica¬ 
tions  are  that  the  two-dimensional  signature  for  the  explosion  source  depends  on 
medium,  and  structure  in  a  complicated  manner.  These  calculations  demonstrate  that 
the  point  source  explosion  model  may  be  inadequate  to  model  discrimination  problems 
in  the  regional  distance  range.  Higher  order  models  for  the  explosion  source  are 
necessary  to  understand  the  effects  of  depth  of  burial,  medium,  and  structure,  on  the 
high/low  frequency  discriminant. 

Recommendations  are  made  to  improve  modeling  of  two-dimensional  explosion 
source  effects  in  short  period  and  high  frequency  regional  waveforms.  Improved  non¬ 
linear  near-source  structure  interaction  and  additional  calculations  for  larger  ranges  of 
scaled  depth  of  burial  and  yield  will  be  required  to  further  test  the  hypothesis  that  the 
empirical  high/low  frequency  regional  discriminant  is  related  to  scaled  depth  of  burial 
effects. 
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II.  INTRODUCnON 

Regional  stations  or  arrays  are  required  for  the  reliable  detection  and  discrimina¬ 
tion  of  small  seismic  events.  At  regional  distances,  small  events  are  recorded  that  do 
not  excite  sufficient  long-period  surface  waves  that  the  traditional  discriminant 

can  be  used.  To  supplant  the  teleseismic  discriminant,  a  variety  of  regional 

discriminants  have  been  suggested  over  the  years  as  reviewed  by  Pomeroy,  et  al. 
(1982)  and  Blandford  (1982).  Tests  of  some  of  these  proposed  regional  discriminants 
have  been  made  in  recent  years  with  western  United  States  earthquakes  and  explo¬ 
sions.  The  most  successful  of  the  regional  discriminants  have  been  regional  variants  on 
the  discriminant  and  a  high/low  frequency  discriminant  (Murphy  and  Bennett, 

1982;  and  recently  Taylor,  et  al.,  1988).  The  empirical  results  of  these  high/low  fre¬ 
quency  discriminants  are  directly  the  opposite  of  the  predictions  of  Evemden,  et  al. 
(1986).  Evemden,  et  al.  suggested  that  explosions  at  regional  distances  will  appear 
enriched  in  "high"  frequencies  relative  to  earthquakes.  The  results  from  careful  studies 
using  Nevada  Test  Site  (NTS)  explosions  and  earthquakes  both  in  the  vicinity  of  NTS 
(Murphy  and  Bennett  1982)  and  earthquakes  dispersed  across  the  western  United 
States  (Taylor,  et  al.  1988)  have  shown  the  opposite  to  be  true.  Earthquakes  have 
appeared  enriched  in  high  frequencies  relative  to  explosions  for  Pn,  Pg,  and  Lg  phases. 
Because  there  is  no  evidence  that  there  are  fundamental  differences  in  the  propagation 
of  the  earthquake  and  explosion  Pn,  Pg,  and  Lg  phases,  any  observed  differences  must 
be  due  to  excitation  differences  between  earthquakes  and  explosions. 

Either  these  differences  reflect  fundamental  spectral  differences  in  the  earthquake 
and  explosion  source,  or  are  due  to  depth  and/or  excitation  differences  in  the  Green’s 
functions.  A  possibly  important  clue  is  the  observation  that  a  few  overburied  explo¬ 
sions  at  NTS  have  high  frequency  enhancement  (Taylor,  et  al.  1988)  and  fail  the 
high/low  frequency  discriminant.  This  suggests  that  scaled  depth  of  burial  effects  par¬ 
tially  control  the  frequency  content  of  regional  phase  excitation.  Taylor,  et  al.  suggest 
that  spall  enhances  the  low  frequencies  of  normal  buried  explosions  and  that  over- 
ouried  explosions  exhibit  less  spall  and  therefore  appear  enhanced  in  high  frequencies. 
Other  possible  explanations  are  (1)  strong  attenuation  of  high  frequencies  from  shallow 
explosions,  or  (2)  differences  between  a  point  explosion  source  model  and  a  more  real¬ 
istic  (2D  or  3D)  model  that  incorporates  asymmetry  in  the  explosion  source. 

In  this  work  we  explore  nonspherically  symmetric  nonlinear  explosion  source 
models  as  functions  of  near  source  structure,  burial  depth,  and  prestress.  Nonlinear, 
two-dimensional,  axisymmetric  finite  difference  calculations  are  studied  to  evaluate  the 
excitation  of  regional  seismic  waves  from  these  calculations.  We  use  calculations  that 
serve  as  explosion  models  in  tuff  and  granite  for  fixed  yield  and  different  depths.  In 
this  way  we  examine  the  relative  frequency  content  of  waves  excited  by  non-point 
source  models  for  explosions.  These  axisymmetric  simulations  represent  the  next 
highest  level  of  complexity  above  the  spherically  symmetric  explosion  point  source. 
Instead  of  simple  linear  "pP"  and  "pS"  reflections  near  the  source,  the  full  nonlinear 
free-surface  interaction  is  simulated  including  nonlinear  yielding,  cracking,  and  spall. 
The  material  properties  responsible  for  the  containment  of  the  explosion  are  functions 
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of  hydrostatic  pressure  and  therefore,  there  is  inherent  symmetry  about  the  vertical  axis 
for  the  simplest  models  that  include  depth  dependent  effects. 

In  addition  to  a  study  of  the  depth  of  burial  effects,  we  study  the  excitation  of  the 
entire  regional  waveform  by  explosions  with  particular  attention  to  the  excitation  of  the 
SV  component  of  motion.  Bennett,  et  al.  (1987)  studied  variations  in  excitation  of 
regional  phases  for  point  explosion  and  earthquake  sources  as  functions  of  depth. 
Although  they  found  that  there  are  significant  variations  in  the  excitation  of  these 
phases  as  functions  of  depth  and  focal  mechanism,  they  found  no  physical  mechanisms 
to  explain  the  empirical  high/low  frequency  discriminant  between  explosions  and  earth¬ 
quakes.  Furthermore,  they  found  as  others  have  noted  that  point  explosion  sources  are 
poor  sources  of  Lg.  Bennett,  et  al.  suggested  that  P-to-SV  conversion  and  "scatter¬ 
ing"  at  interfaces  could  produce  the  required  SV  energy  to  couple  into  the  Lg 
waveguide.  It  is  well  known  that  explosion  generated  Lg  is  a  robust  measure  of  the 
explosion  size  (Nuttli,  1986)  in  contrast  to  this  theoretical  difficulty  in  modeling  the 
phase.  Patton  (1988)  in  a  detailed  study  of  the  HARZER  explosion  suggests  that  spall 
is  a  significant  contributor  to  regional  Lg.  Using  the  nonlinear  2D  finite  difference 
simulations  we  test  the  hypothesis  that  axisymmetric  explosion  models  produce  more 
SV  than  point  explosion  sources. 

We  examine  the  differences  between  the  point  explosion  source,  and  explosion 
models  derived  from  nonlinear,  two-dimensional,  axisymmetric  finite  difference  calcu¬ 
lations.  Although  these  calculations  do  not  simulate  azimuthal  variability  in  the 
source,  they  do  simulate  nonlinear  free  surface  interactions  including  spall,  and  asym¬ 
metries  that  depend  on  depth  of  burial.  These  models  serve  as  a  "first  principles" 
approach  that  are  a  step  more  complicated  than  spherically  symmetric  explosions.  In 
the  case  of  cratering  simulations,  even  the  physics  of  crater  formation  arc  modeled. 

The  work  reported  here  is  divided  into  three  main  sections.  These  sections,  (1) 
introduce  a  velocity  model  for  the  propagation  of  regional  waves  in  the  basin  and 
range,  (2)  examine  far-field  P  and  SV  waves  radiated  by  axisymmetric  nonlinear  finite 
difference  simulations  of  explosions,  and  (3)  introduce  a  new  method  for  the  propaga¬ 
tion  of  the  nonlinear  explosion  simulations  to  regional  distances.  A  brief  discussion  of 
the  conclusions  and  recommendations  follows. 

In  Section  III,  we  present  a  velocity  model  for  purposes  of  propagation  of 
regional  seismic  signals,  and  examine  some  results  of  point  source  synthetics  of  explo¬ 
sions  and  earthquakes.  These  calculations  serve  two  purposes.  First,  we  examine  a 
velocity  and  Q  model  for  the  Basin  and  Range  based  on  a  literature  survey  and  vali¬ 
date  the  model  by  examining  the  distance  dependence  of  Pn,  Pg,  and  Lg.  Second,  we 
confirm  the  conclusions  of  Bennett,  et  al.  (1987)  that  the  regional  propagation  and 
excitation  differences  between  earthquakes  and  explosions  does  not  alone  explain  the 
high/low  frequency  discriminant.  As  part  of  a  validation  and  testing  process  we  have 
examined  regional  short-period  and  broad-band  seismograms  in  the  Western  United 
States  in  order  to  compare  their  propagation  character  (spectral  and  time  domain)  to 
predictions  from  Green’s  functions  produced  by  the  wavenumber  integration  code.  We 
review  results  of  several  authors  for  the  attenuation  of  regional  phases  in  the  western 
United  States.  Procedures  similar  to  those  applied  to  real  data  are  then  applied  to 
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synthetic  Green’s  functions  to  validate  the  velocity  and  attenuation  models  for  western 
United  States  regional  wave  propagation.  Results  indicate  that  frequency  dependent  Q 
models  may  be  required  for  the  crust  in  order  to  adequately  explain  the  frequency 
dependence  of  Pg  and  Lg  as  a  function  of  distance.  This  warrants  additional  study  for 
the  understanding  of  regional  discriminants. 

In  Section  FV,  we  present  far-field  P  and  SV  waveforms  from  axisymmetric  non¬ 
linear  finite  difference  calculations  for  Pahute  Mesa  and  Shagan  River  structures  at 
various  scaled  depths  of  burial.  In  addition,  far-field  P  and  SV  from  nonlinear  models 
for  PILEDRIVER  with  and  without  prestress  are  compared.  We  have  extended  the 
work  of  Day,  et  al.  (1986b)  that  compared  predicted  teleseismic  P-waveform  statistics 
to  observed  statistics  of  Log(Pmax/Pa)  and  Log(Pb/Pmax)  of  McLaughlin,  et  al. 
(1985).  For  comparison  with  results  of  Der,  et  al.  (1987),  McLaughlin,  et  al.  (1986) 
and  McLaughlin,  et  al.  (1987)  we  have  generated  a  suite  of  synthetics  based  on  a 
number  of  nonlinear  calculations  intended  to  simulate  Shagan  River  and  Pahute  Mesa 
contained  and  cratering  explosions  as  well  as  contained  explosions  in  granite.  These 
synthetics  are  compared  to  observed  seismogram  statistics.  Using  adjustments  for 
depth  of  burial  based  on  synthetic  waveforms,  the  150  kt  limit  is  estimated  for  the 
Shagan  River  test  site. 

We  conclude,  that  in  general,  the  2D  nonlinear  calculations  under  predict  the 
short  period  teleseismic  "pP"  amplitude.  The  2D  nonlinear  calculations  either  over 
predict  the  amount  of  attenuation  of  the"pP’'  phase  or  under  predict  the  nonlinear  spall 
contribution  to  this  phase.  The  P-wave  spectral  scaling  of  the  2D  nonlinear  calcula¬ 
tions  suggest  that  overburial  enhances  high  frequencies  with  respect  to  normal  burial. 
However,  2D  nonlinear  SV  radiation  is  enhanced  in  low  frequencies  relative  to  that 
predicted  by  a  point  source  model.  Taken  together,  the  P-wave  and  SV-wave  radiation 
predicted  by  the  2D  nonlinear  calculations  are  inconsistent  with  a  simplified  ID  point 
source  model  for  the  explosion. 

Finally  in  Section  V,  we  propagate  the  axisymmetric  nonlinear  simulations  to 
regional  distances  using  a  newly  developed  interface  between  the  finite  difference 
codes  and  a  wavenumber  integration  synthetic  code.  The  theoretical  justification  based 
on  the  elastodynamic  representation  theorem  for  this  interface  is  contained  in  an 
appendix.  Finite  difference  explosion  simulations  are  propagated  to  regional  distances 
by  the  use  of  the  wavenumber  integration  code.  This  simulates  aspects  of  the  source 
finiteness  and  directivity  not  modeled  by  explosive  point  sources.  In  this  manner,  the 
two-dimensional  aspects  of  the  source  models  can  be  properly  propagated  to  regional 
distances. 

We  find  that  enhanced  shear  wave  (Lg)  generation  is  evident  in  some  axisym¬ 
metric  simulations  due  to  the  nonlinear  free  surface  interaction.  This  addition  to  the 
point  explosion  source  model  can  to  first  order  be  seen  as  a  vertically  oriented  com¬ 
pensated  linear  vector  dipole  (CLVD)  source. 

Although  the  2D  nonlinear  synthetics  qualitatively  show  more  explosion  gen¬ 
erated  Lg  than  the  ID  point  source  model,  the  Lg  excitation  is  still  well  below  that 
commonly  observed  from  explosions.  The  2D  nonlinear  synthetics  suggest  that 
processes  that  enhance  the  CLVD  component  of  the  explosion  source  will  contribute  to 
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regional  SV  energy.  Such  processes  include,  but  are  not  limited  to,  spall  and  near 
source  scattering.  Simuladons  of  explosions  in  prestressed  granite  show  remarkably 
little  difference  between  no  prestress  simulations  in  the  short-period  frequency  range. 
These  conclusions  are  in  agreement  with  the  lack  of  correlation  of  short  period  (Lg) 
observations  with  long  period  (LR)  tectonic  release  observations. 

Day,  et  al.  (in  preparation)  compare  two  popular  methods  for  the  construction  of 
synthetic  regional  seismograms.  The  method  of  modal  summation  with  "locked- 
modes"  as  described  by  Harvey  (1981)  is  compared  to  a  wavenumber  integration  pro¬ 
gram  based  on  the  work  of  Apsel  (1979).  It  is  found  that  the  two  methods  produce 
consistent  results  for  high  Q  Earth  models.  However,  the  perturbation  approximation 
inherent  in  the  modal  summadon  procedure  can  break  down  for  low  Q’s  appropriate 
for  the  western  United  States.  Therefore,  the  wavenumber  integradon  method  for  com- 
putadon  of  synthedc  regional  seismograms  is  the  most  accurate  method  if  attenuadon 
is  to  be  considered.  Also,  modal  summadon  computadon  time  rapidly  increases  with 
increasing  frequency,  such  that  the  wavenumber  integradon  procedure  is  faster  and 
more  accurate  at  frequencies  above  1  Hz.  Therefore,  except  where  stated,  the 
wavenumber  integradon  method  is  used  for  the  computadon  of  regional  synthedcs 
throughout  this  report. 
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in.  POINT  SOURCE  REGIONAL  SYNTHETICS 

3.1  Introduction 

In  order  to  properly  model  regional  seismograms,  we  require  a  model  for  the 
crust  and  upper  mantle  with  appropriate  velocity  and  attenuation  properties.  In  this 
section  we  introduce  a  velocity  and  Q  model  for  the  basin  and  range  province  and 
examine  the  behavior  of  explosion  and  earthquake  Green’s  functions  for  this  model. 
The  model  is  chosen  with  a  frequency  independent  Q  in  crust  and  upper  mantle  and  is 
validated  by  comparison  of  the  spatial  decay  of  regional  phases  for  synthetics  and 
observations.  These  calculations  provide  a  baseline  for  the  evaluation  of  regional  syn¬ 
thetics  based  on  the  nonlinear  two-dimensional  finite  difference  source  simulations. 

3.2  Basin  and  Range  Model,  BR2 

The  velocity  and  Q  model,  BR2  listed  in  Table  3.1  (Figures  3.1a,b),  was  intended 
to  be  representative  of  a  number  of  velocity  models  that  have  been  proposed  for  the 
basin  and  range  province  of  the  western  United  States.  Slight  gradients  were  intro¬ 
duced  in  the  upper  crust,  lower  crust  and  upper  mantle  lid  in  order  to  turn  rays  but  not 
produce  noticeable  triplications.  An  upper  mantle  lid  P-wave  velocity  gradient  was 
chosen  to  be  just  below  that  of  the  CITKX)  model  series,  and  the  Herrin,  et  al.  (1968) 
P-wave  models  (0.01215  km/s/km).  The  velocity  model  is  not  intended  to  represent 
any  single  path  but  to  be  representative  of  a  typical  Basin  and  Range  path.  The 
CITKX)  series  models  of  Archambeau,  Flinn,  and  Lambert  (1969),  as  well  as  the 
models  of  Pakiser,  (1963),  Helmberger  (1973),  Priestley  and  Brune  (1978),  Priestley, 
et  al.  (1982),  Patton  and  Vergino  (1981),  Chun  (1983),  Wallace  (1983),  Patton  and 
Taylor  (1984),  Taylor  and  Patton  (1986),  and  the  NTS4  model  of  Bennett,  et  al. 
(1987)  were  all  considered.  It  was  decided  to  suppress  the  crustal  low- velocity  zones 
of  the  models  and  retain  only  the  well  established  low  velocity  zone  of  the  upper  man¬ 
tle.  There  is  a  strong  initial  velocity  gradient  at  the  top  of  the  crust  and  a  lesser  gra¬ 
dient  throughout  the  rest  of  the  crust. 
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Figure  3.1b.  Basin  and  Range  Q  model,  BR2. 
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TABLE  3.1  Model  BR2 


h 

(km) 

a 

(km/s) 

(km/s) 

P 

(gm/cc) 

Qp 

Qs 

0.50 

3.00 

1.58 

2.00 

150. 

75. 

0.50 

3.80 

2.20 

2.20 

200. 

100. 

2.00 

5.80 

3.40 

2.80 

250. 

125. 

3.00 

6.00 

3.46 

2.90 

250. 

125. 

14.0 

6.30 

3.64 

2.90 

300. 

150. 

4.00 

6.50 

3.75 

2.90 

300. 

150. 

4.00 

6.80 

3.93 

3.00 

300. 

150. 

6.00 

6.90 

3.98 

3.00 

300. 

150. 

4.00 

7.80 

4.50 

3.20 

400. 

200. 

4.00 

7.84 

4.50 

3.20 

400. 

200. 

4.00 

7.88 

4.50 

3.20 

400. 

200. 

4.00 

7.92 

4.50 

3.20 

400. 

200. 

4.00 

7.96 

4.50 

3.20 

400. 

200. 

2.00 

8.00 

4.50 

3.20 

400. 

200. 

4.00 

8.05 

4.50 

3.20 

400. 

200. 

10.0 

7.90 

4.40 

3.20 

150. 

75. 

10.0 

7.80 

4.20 

3.20 

150. 

75. 

20.0 

7.80 

4.20 

3.20 

150. 

75. 

20.0 

7.90 

4.50 

3.20 

150. 

75. 

60.0 

8.00 

4.55 

3.30 

150. 

75. 

20.0 

8.10 

4.60 

3.40 

300. 

150. 

The  Q  model  also  was  intended  to  be  representative  of  the  basin  and  range  but 
not  to  represent  any  one  particular  path.  Although  there  may  be  some  recent  evidence 

Qp 

for  deviation  from  the  rule  that  =  2,  the  model  was  constructed  with  this  con- 

Qs 

straint.  Crustal  shear- wave  estimates  vary  quite  significantly  from  author  to  author 
and  much  of  the  variation  may  simply  be  lateral  variation  in  attenuation.  Also,  fre¬ 
quency  dependence  of  Q  is  evident  in  the  data  and  produces  differences  between 
authors.  Consequently,  a  median  shear-wave  average  Q  value  near  1  Hz  was  chosen 
for  the  crust  and  several  depth  dependent  Q  models  were  used  for  guidance  to  choose 
a  depth  dependence.  Depth  dependent  Q  models  for  the  region  include  Bache,  et  al. 
(1980),  Chun  (1983),  Patton  and  Taylor  (1984),  and  Taylor  and  Patton  (1986).  The 
chosen  model  features  a  low  surface  Q  in  order  to  suppress  short-period  fundamental 
mode  at  large  regional  distances  and  a  low  Q  zone  in  the  upper  mantle  low-velocity 
zone.  Q  estimates  for  Pn  indicate  that  there  exists  a  high  Q  lid  to  the  mantle  low- 
velocity  zone.  A  suimnary  of  crustal  and  Qig  estimates  from  the  literature  is  tabu¬ 
lated  below  (Table  3.2)  and  graphically  portrayed  in  Figure  3.2. 

In  general  however,  the  BR2  model  may  actually  underestimate  the  total  amount 
of  attenuation  in  the  crust  and  upper  mantle  at  1  Hz.  Model  BR2  has  a  teleseismic  P- 
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Figure  3.2. 


Crustal  Q-  and  estimates  from  the  literature  as  a 

function  m  frequency.  Lg  studies  are  designated  with 
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wave  t*  of  about  0.15  second,  and  a  teleseismic  S-wave  t^of  0.50  second  down  to  a 
depth  200  km.  This  is  0.20  seconds  less  than  the  total  t*p  that  Der,  et  al.  (1985) 
would  attribute  to  the  NTS  source  region,  and  would  leave  the  remainder  to  be 
accounted  for  by  attenuation  along  the  remainder  of  the  mantle  P-wave  path  to  the 
turning  point  of  the  ray. 

The  average  Q’s  of  the  BR2  model  cmst  averaged  over  vertical  travel  time  for  P 
and  S  waves  are  139  and  279,  respectively.  The  average  crustal  S-wave  Q  of  about 
140  was  chosen  to  coincide  with  the  median  value  for  Lg  Q’s  at  1  Hz  from  a  number 
of  researchers  working  in  the  Basin  and  Range. 


TABLE  3.2  SUMMARY  OF  AVERAGE  CRUSTAL  Qr  ESTIMATES 
_ REFERENCE _ Q _ RANGE  (Hz) 


Press  (1964) 

Mitchell  (1975) 

Lee  and  Solomon  (1975) 

Braile  (1977) 

Bache  et  al.  (1978)(NTS-ALQ) 
Cheng  and  Mitchell  (1981) 

Deret  al.  (1981)(OB2NV) 

Chun  (1983) 

Singh  and  Herrmann  (1983) 

Peseckis  and  Pomeroy  (1984) 

Patton  and  Taylor  (1984) 

Nuttli  (1981) 

Nuttli  (1986)(NTS-TUC-DUG-BKS) 
Chavez  and  Priestley  (1986)(EXP) 
Chavez  and  Priestley  (1986)(EQKS) 
Taylor  et  al.  (1988) 

Patton  (1987)(NTS-ELK) 

Patton  (1987)(NTS-KNB) 

Patton  (1987)(NTS-LAC) 

Patton  (1987)(NTS-MNV) 


Qig 

=  450 

0.7- 1.3 

Qp  = 

125 

0.05  to  0.1 

(2p  = 

200 

0.05  to  0.1 

60 

=  1 

(2p  = 

=  125 

0.05  to  0.1 

Gp  = 

85 

0.2 

Qig  '■ 

=  200 

0.5  to  8 

42 

0.05  to  0.2 

QLg 

=  250/°'‘5 

0.5  to  2 

Qig 

=  264/0-^ 

0.5  to  2 

ep  = 

110 

0.02  to  0.15 

Qig 

=  290/®'^ 

1  to  5 

QLg 

=  290/®"^ 

0.3  to  2 

QLg 

=  206/®-^* 

0.3  to  10 

QLg 

^214/®-54 

0.3  to  5 

QLg 

=  150/®-^ 

1  to  6 

QLg 

=  150/®-^ 

0.3  to  3.0 

QLg 

=  142/®-^ 

0.3  to  3.0 

QLg 

=  97/®'^ 

0.3  to  3.0 

QLg 

=  93/®-^ 

0.3  to  3.0 

The  upper  mantle  Q  model  was  chosen  to  strike  a  balance  between  relatively  high 
Q’s  observed  for  Pn  and  the  surface  wave  models  that  have  less  resolution  as  to  the 
location  of  the  low-Q  region  in  the  upper  mantle.  Estimates  of  the  average  Q  for  Pn 
are  often  dependent  on  the  geometrical  spreading  assumed  and  hence  are  velocity 
model  dependent.  A  model  that  assumes  a  half-space  upper  mantle  will  predict  a  fre¬ 
quency  dependent  Pn  propagation  with  the  classical  head-wave  nature.  Observations 
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of  Pn,  however,  show  that  it  is  an  impulsive  high-frequency  arrival  and  that  the  Pn 
spectrum  corrected  for  the  source  spectrum  does  not  fall-off  fester  than  proportional  to 

r\ 

In  order  to  demonstrate  feis,  we  correct  explosion  Pn  spectra  for  instr^ent  and 
source  and  then  estimate  a  t*pn  from  the  spectral  slope.  Estimates  of  for  a 
number  of  paths  are  listed  below  (Table  3-3)  for  the  HARDHAT  explosion  (Climax 
Stocl^NTS)  to  several  regional  LRSM  (Long  Range  Seismic  Measurement)  stations. 
The  estimates  would  indicate  a  mean  Q  of  about  6(X)  with  a  standard  error  of  the 
mean  of  100  and  a  standard  error  of  the  observations  of  400.  This  large  scatter 
implies  that  individual  event-receiver  paths  exhibit  large  variations  in  attenuation. 
Apart  from  the  large  scatter  in  these  estimates  the  mean  value  is  in  general  agreement 
with  an  apparent  (2  a  around  500  in  the  upper  mantle  lid  if  we  account  for  the  average 
attenuation  of  the  crust. 


TABLE  3.3 

HARDHAT  TO  LRSM 

(0.5  to  7  Hz) 

A(deg) 

t*P„(sec) 

STATION 

1.4 

0.06 

DV-CL 

2.6 

0.05 

KN-UT 

3.0 

0.05 

TN-CL 

3.6 

0.14 

WM-AZ 

3.6 

0.05 

FM-UT 

4.3 

0.05 

'VI-NV 

4.4 

0.14 

FS-AZ 

4.5 

0.40 

CP-CL 

5.3 

0.05 

SF-AZ 

6.6 

0.10 

DR-CO 

6.6 

0.20 

HL-ID 

8.7 

0.40 

PT-OR 

9.2 

0.85 

LC-NM 

9.4 

0.40 

RT-NM 

All  of  the  CmOO  series  models  of  Aichambeau,  Flinn,  and  Lambert  (1969)  had 
very  high  Q  lids.  Werth,  Herbst,  and  Springer  (1962)  estimated  =  300  to  400  in 
the  96  to  714  km  distance  range  based  on  observations  of  Ranier  Mesa  events  recorded 
on  Benioff  seismograph  systems.  At  the  high  end  of  estimates  from  the  literature, 
Evemden,  et  al.  (1986)  suggest  that  Pn  has  a  Q  between  1000  and  2000  in  tectonic 
regions. 

As  far  as  frequency  dependence  for  model  Q’s  are  concerned,  we  have  chosen  to 
make  the  synthetic  computations  (at  least  at  first)  with  a  frequency  independent 
attenuation  structure.  However,  numerous  authors  have  reported  frequency  dependent 
Q  models  for  Lg  as  listed  above,  fewer  authors  have  reported  frequency  dependence 
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for  Pn  and  Pg  attenuation.  Taylor,  et  al.  (1988)  estimate  frequency  dependence  for 
Pn  and  Pg  apparent  Q’s  as  Qp^  =  250  and  Qp.  =  300  ^  These  correspond  to 

—  *  Ht* 

a  1  Hz  i*p^  at  300  km  of  0.15  sec  and  a  r*p„  =  t*p„  +  f  — —  =  0.06  sec.  The  effect 

—  4/ 

of  the  power  law  frequency  dependence  is  to  lower  t*p„  which  controls  the  spectral 
shape  of  the  phase  and  if  interpreted  as  a  constant  Q  model  would  imply  a  Qp^  in 
excess  of  600.  Similarly  at  300  km  and  1  Hz,  Taylor,  et  al.  would  predict  t*pg  = 
0.14  sec. 


Qa 

We  have  assumed  that  =  2.  Both  analysis  presented  later  in  this  section  and 

Q  p 

the  results  of  the  Q  estimates  of  Taylor,  et  al.  (1988)  suggest  mat  for  frequencies 
above  1  Hz,  the  crustal  P-  and  S-wave  Q’s  may  have  different  frequency  dependence 
and  that  at  high  frequencies,  may  be  higher  than  Similar  results  were  found 
by  Goncz,  et  al.  (1986)  for  the  eastern  United  States  for  a  comparison  of  apparent  Pn 
and  Sn  Q’s.  Also,  Soreno  and  Orcutt  (1987)  claim  similar  results  for  oceanic  Pn  and 
Sn  phases.  It  is  not  known  at  this  time  if  this  is  due  to  propagation  effects  in  layered 
media,  intrinsic  bulk-attenuation,  or  whether  this  is  due  to  scattering  by  lateral  varia¬ 
tions  in  crastal  structure  producing  different  effects  on  P  type  and  S  type  modes  of 
propagation.  From  our  numerical  experiments,  it  appears  that  this  phenomenon  is  not 
due  to  one-dimensional  (ID)  layered  structure  effects,  which  leaves  intrinsic  bulk- 
attenuation  and/or  scattering.  Extensions  to  the  wavenumber  integration  procedure  are 
easily  implemented  that  will  introduce  separate  frequency  dependence  into  the  attenua¬ 
tion  of  P  and  S  waves  in  the  medium. 


3.3  Amplitude-Distance  Relationships 

We  also  wish  to  represent  the  average  observed  amplitude-distance  relationships 
for  Pn,  Pg,  and  Lg.  In  the  following,  we  review  some  of  the  sometimes  conflicting 
reports  from  the  literature.  These  differences  often  reflect  the  difference  between  using 
Log(amplitude/period),  versus  Log(amplitude),  or  maximum  amplitude  versus  a 
specified  cycle,  or  differences  in  instrumental  bandwidths. 

As  far  as  empirical  distance  relationships  are  concerned,  the  Pn  amplitudes  from 
explosions  or  shallow  earthquakes  in  the  western  U.S.  are  generally  proportional  to 
to  according  to  a  review  by  Blandford,  et  al.  (1981).  The  Werth,  et  al. 
(1962)  data  can  be  regressed  to  a  relationship  for  three  events  (10  data  points)  as 
Log(Pn)  =  Log(Yield)  +  9.3(0.6)  -  3.05(0.24)*Log(R),  for  the  amplitude  of  the  first 
half  cycle  of  the  Pn  (amplitude  in  nm,  R  in  km).  Cases  for  attenuation  less  than  R~^ 
may  be  found,  such  as  the  Veith  and  Clawsen  (1972)  B  values  for 
Log(amplitude/period)  zero  to  peak  P-wave  magnitude  which  correspond  to  a  linear 
formula  for  mf,  from  A  =  2  to  10  degrees, 

rrih  =  -1.55(0.08)  +  2.1(0.1)Log  (A). 

Evemden  (1967)  derived  a  formula  for  from  Pn  peak  to  peak  amplitude  (in  nm)  of 
mi,  =  -7.55  -I-  \.2lLog{A/T)  +  2>.QALog{R). 
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The  Evemden  study  was  based  on  LRSM  records  for  NTS  explosions.  Taylor,  Denny, 
and  Vergino  (1986)  found  using  broadband  velocity  seismograms  that  Lg(A/T)  was 
proportional  to  R~^  for  earthquakes  and  explosions  in  the  300  to  1000  km  range. 
However  Denny,  Taylor  and  Vergino  (1987)  arrived  at  a  preferred  formula  of 

=  -3.95  +  Log  (A )  +  lAlLog  {R ) 

and  a  distance  coefficient  of  2.79  for  Log(A/T).  They  argued  that  mi,  based  on 
Log(A)  was  more  stable  than  that  based  on  Log(A/T)  and  that  better  :Ms  discrimi¬ 
nation  of  earthquakes  and  explosions  resulted  from  their  choice  for  the  mi,  distance 
correction  formula.  Clearly,  we  see  that  instrumental  bandwidths,  and  the  choice  of 
Log(A)  versus  Lx)g(A/r)  are  important  factors  in  the  published  amplitude  distance 
decays  for  Pn. 

These  differences  in  instrumental  bandwidth  impact  the  use  of  other  phases  as 
well  as  Pn.  Blandford,  et  al.  (1981)  report  that  Pg  or  Pmax  as  measured  on  LRSM 
stations  decays  as  R~^^  in  the  western  U.S.  Similarly,  they  found  that  Lg  decayed  as 
Based  on  these  distance  corrections  they  reported  that  the  distance  corrected 
Pg/Lg  ratio  could  be  used  to  discriminate  between  earthquakes  and  explosions.  Barker 
(1981)  using  raw  ANMO  records  of  NTS  explosions  and  various  WUS  earthquakes 
also  reported  separation  of  earthquake  and  explosion  populations  using  the  Pg/Lg  ratio. 
Barker’s  explosions  were  of  varying  size  but  at  nearly  the  same  distance  range. 
Barker  also  reponed  success  using  Pg/Pn  on  the  raw  records  and 

Pn(0.625Hz)/Pn(2,5Hz)  ratios  in  conjunction  with  the  Pg/Lg  ratio.  Murphy  and  Ben¬ 
nett  (1982)  using  earthquakes  and  explosions  in  southern  Nevada  recorded  at  TFO 
found  the  Pg/Lg  ratios  from  the  raw  records  did  not  discriminate  eanhquake  and 
explosion  populations.  However,  they  found  that  Lg(2-4Hz)/Lg(0.5-lHz)  ratios  suc¬ 
cessfully  separated  explosions  and  earthquakes.  Gupta,  et  al.  (1984)  reported  that 
high-fi^uency  to  low-frequency  ratios  for  Lg  were  a  useful  discriminant,  but  that  a 
spectral  measure  of  Pg/Lg  could  also  be  used  to  increase  the  discrimination  power  of 
the  two  combined  statistics.  Since  these  studies  were  conducted  with  different  instru¬ 
mental  bandwidths  (LRSM,  SRO,  Johnson-Matheson,  and  Halls-Sears),  the  frequency 
and  distance  dependence  of  raw  record  Pg/Lg  ratios  may  have  contributed  to 
discrepancies  between  some  of  the  workers.  The  combination  of  two  or  more  meas¬ 
ures  of  the  regional  phases  in  a  discriminant  function  may  overcome  some  of  the  vari¬ 
ations  in  a  single  discriminant  measure.  Clear'),  if  these  difficulties  are  to  be  over¬ 
come,  then  these  phase  ratio  and  spectral  disc’^  ninants  must  take  into  account  the  dis¬ 
tance  dependence  of  the  phases  in  order  to  make  frequency  and  distance  corrections. 

Taylor,  et  al.  (1988)  extended  the  'lurphy  and  Bennett  discriminant  to  higher 
frequencies  and  found  that  the  (6-8Hz)/(l-2Hz)  works  even  better  than  the  (2- 
4Hz)/(0.5-lHz)  measure.  In  order  to  accomplish  this,  they  introduced  distance  and  fre¬ 
quency  dependent  corrections  for  each  phase.  However,  they  did  find  exceptions  for 
some  overburied  explosions.  The  overburied  explosions  appear  to  have  higher  fre¬ 
quency  Lg  than  normally  contained  explosions. 

In  conclusion,  we  see  that  in  order  to  make  a  portable  discriminant  based  on  the 
spectra  of  single  phases  or  combinations  of  phases,  will  require  an  understanding  of 
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the  basic  propagation  characteristics  for  the  regional  wavepackets.  In  the  next  section, 
we  show  that  the  full  waveform  synthetics  can  simulate  these  characteristics  and  then 
discuss  some  characteristics  of  the  model  required  to  achieve  this. 

3.4  Analysis  of  BR2  Synthetics 

Figtu^s  3.3a,  b,  and  c  show  synthetics  at  325  km  and  400  km  distance  compared 
to  records  from  the  broadband  stations  KNB  (*320),  LAC  ("320),  and  ELK  (~400)  for 
BACKBEACH  and  MAST.  The  synthetics  are  velocity  in  nm/sec  for  a  Heavyside 
step  explosion  point  source  (moment  of  10*^  Nt-M)  at  a  depth  of  1  km  for  the  BR2 
mt^l.  The  data  is  recorded  on  the  LLNL  broadband  systems  with  a  flat  response  to 
velocity  in  the  bandwidth  of  interest.  The  data  and  synthetics  have  been  bandpassed  in 
four  bands  0.05-0.16  Hz,  0.16-0.50  Hz,  0.50-1.58  Hz,  and  1.58-5.0Hz.  The  general 
character  of  the  Pn  and  Pg  phases  at  high  frequencies  and  the  low-frequency  surface 
waves  are  similar  to  the  real  data.  The  primary  differences  can  be  seen  in  the  coda  of 
the  deterministic  wavepackets  and  in  the  amplitudes  of  the  Lg  at  intermediate  frequen¬ 
cies  and  above  1  Hz.  BACKBEACH  can  be  seen  to  have  generated  relatively  more 
Lg  than  MAST  at  ELK  and  KNB.  However,  in  general  the  Lg  amplitudes  are  not 
more  than  a  factor  of  two  above  the  Pg  coda.  Qearly,  the  explosion  point  source  does 
not  excite  enough  Lg  relative  to  other  phases  to  agree  with  observations.  Also,  the 
relative  amplitude  of  the  synthetic  Lg  signals  between  the  0.50-L58Hz  and  the  1.58- 
5.0Hz  bandwidths  is  not  right.  The  explosion  Lg  that  is  excited  in  the  ID  layered 
models  by  a  point  source  does  not  contain  enough  energy  above  1  Hz.  We  will  return 
to  this  later  after  a  more  detailed  analysis  of  the  synthetics  using  both  spectral  and 
time  domain  bandpass  analysis.. 

In  order  to  test  that  the  velocity  and  Q  model  represents  propagation  in  the 
western  U.S.  we  have  used  the  synthetics  as  data  is  often  used  to  infer  the  spatial 
decay  of  regional  phases  and  estimate  Q.  Using  broadband  velocity  records  that 
represent  the  velocity  response  to  a  step  function  source  we  find  that  for  a  deviatoric 
source  (average  between  strike-slip  and  dip-slip  at  45°  azimuth)  at  8  km  depth  the 
peak  amplitudes  (velocity)  of  the  three  phases  can  be  characterized  as: 

Log(Pnmax)  «  -2.84*Log(R/l()0),  o  =  0.14 

Log(Pgmax)  «  -2.40*Log(R/100),  a  =  0.13 

Log(Lgmax)  «  -2.44*Log(R/100),  a  =  0.10 

The  power  law  decay  of  these  phases  between  R~^  and  are  in  rough  agreement 
with  other  broadband  velocity  decays  observed  in  the  region.  Since  the  ampUtude- 
distance  relationship  is  frequency  dependent,  and  thuefore  will  also  be  dependent  on 
instrumental  bandwidth,  we  applied  two  different  techniques  to  estimate  the  apparent 
spatial  Q  from  the  synthetics. 

3.4.1  Spectral  Domain  Measures  of  Apparent  Q 

We  assume  that  the  amplitude  of  a  phase  is  proportional  to  a  geometrical  spread¬ 
ing  factor  R~^  where  0.5  ^  g  $  1,  and  an  anelastic  factor  expir^fRIclQ)  where  R  is 
distance,  c  is  the  group  velocity,  and  Q  is  the  apparent  quality  factor.  Therefore,  we 
may  estimate  Q  by  either  assuming  a  geometrical  spreading  exponent  and  fit  the 
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Figure  3.3b.  Comparison  of  bandpass  filtered  synthetics  at  325  km  with  observations  at  station 

LAC:  BACKBEACH  (right),  synthetic  (center),  and  MAST  (left). 
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Figure  3.3c.  Comparison  of  bandpass  filtered  synthetics  at  400  km  with  observations  at  station 

ELK;  BACKBEACH  (right),  synthetic  (center),  and  MAST  (left). 
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amplitude  decay  with  distance,  or  we  may  compute  ratios  of  spectra  at  different  dis¬ 
tances  and  estimate  Q  from  the  slope  of  the  spectral  ratio:  Log{A-^IA'^  =  -g 
I^g(/?  xIR  2)  -  nfLog  (e  )iR  y-R  j)lc  IQ . 

Using  spectral  ratios  of  vertical  components  between  different  distances 
(200km/600km,  300km/700km,  etc...)  for  synthetic  Pg  (c  =  6  km/s)  we  find  no  dis¬ 
tance  dependence  for  the  apparent  spectral  Q  and  only  insignificant  differences 
between  the  Pg  excited  by  an  explosion  and  that  by  a  deeper  "earthquake"  (All  "earth¬ 
quake"  results  quoted  here  are  an  average  between  the  response  of  strike-slip  and  dip- 
slip  sources  at  45°  azimuth).  We  find  Qpg  =  256(32)  mean,  270  rms,  231  median,  and 
244(11%)  geometrical  mean  for  the  8  km  deep  earthquake  source.  Similarly  for  the  1 
km  explosion  source,  we  find  Qpg  =  306(50)  mean,  335  rms,  and  276(17%)  geometri¬ 
cal  mean.  The  spread  in  Q  by  either  determination  is  about  +/-  140.  So  even  with 
synthetic  data,  it  is  clear  that  such  a  determination  is  subject  to  considerable  variance. 
For  comparison,  the  BR2  model  has  an  average  Qp  of  about  270  in  the  crust  which 
compares  favorably  with  the  spectral  Q  observed  in  the  synthetic  Pg’s.  We  see  no  dis¬ 
tance  dependence  on  the  apparent  spectral  Q  to  suggest  for  example  that  ratios 
between  200  km  and  400  km  are  any  different  than  ratios  between  4()0  km  and  600 
km. 

Similar  results  for  Pn  apparent  Q  (c  =  7.8  km/s)  from  spectral  ratios  show  a  very 
large  spread  with  Qp^  =  ^7(150)  mean,  743  rms,  370  median,  and  506(22%) 
geometrical  mean.  The  large  uncenainty  from  synthetic  data  indicate  that  the  pros¬ 
pects  of  determining  Q  in  this  manner  from  real  data  will  encounter  even  larger 
scatter.  The  values  are  roughly  consistent  with  the  Qp  value  of  400  in  BR2  for  the 
upper  mande  lid,  but  suggest  that  the  estimates  are  biased  on  the  high  side. 

For  Lg  (c  =  3.5  km/s)  we  find  Q^g  =  242(20)  mean,  250  rms,  253  median  ,and 
240(9%)  geometrical  mean,  for  the  8  Ion  deep  earthquake,  Qig  =  144(13),  150  rms, 
142  median,  140(10%)  geometrical  mean,  for  a  0.5  km  deep  earthquake  source,  and 
Qig  =  106(7)  mean,  108  rms,  102  median,  and  104(4%)  geometrical  mean  for  a  1  km 
deep  earthquake  source.  We  find  no  apparent  distance  dependence  on  the  apparent  Lg 
attenuation.  The  estimates  from  the  shallow  earthquake  synthetics  are  consistent  with 
the  average  Q  for  the  model  crust  of  BR2  of  139.  The  deep  earthquake  Lg  spectra  are 
contaminated  by  Sn  which  arrives  within  the  Lg  window  and  accounts  for  the  higher  Q 
estimate.  The  shallow  earthquake  Lg  is  better  excited  and  the  Q  is  more  consistent 
with  the  average  crustal  Qp. 

3.4.2  Time  Domain  Measures  of  Apparent  Q 

The  time  domain  measures  for  attenuation  are  estimated  by  bandpass  filtering  the 
records  and  then  fitting  a  functional  form  to  the  Log-peak-amplitude  in  each  band  as  a 
function  of  distance.  Trade-off  between  the  geometrical  spreading  assumed  and  the 
apparent  Q  is  unavoidable  in  this  procedure.  However,  bounds  can  be  placed  on  the 
geometrical  effects  that  can  be  expected  for  waveguide  phases  such  as  Pg  and  Lg. 
Geometrical  effects  can  be  bounded  by  cylindrical  spreading,  and  spherical 

spreading,  R~^.  The  Airy  phase  approximation  often  used  for  Lg  results  in  R~^'^ 
(Nuttli,  1986).  Since  regional  attenuation  is  often  described  simply  as  a  power  law, 
we  will  first  simply  give  the  best  fit  exponents  for  each  bandwidA  of  interest.  We 
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have  used  Butterworth  bandpass  filters  of  one  octave  width  and  48  db/octave  roll-offs 
outside  the  bandwidth  specified  by  -3  db  comer  frequencies. 

Using  the  1  km  deep  explosion  and  the  8  km  deep  earthquake  vertical  velocity 
records,  we  can  summarize  the  data  in  the  following  tables.  For  the  three  bandwidths, 
0.625  to  1.25,  1.25  to  2.5  and  2.5  to  5  Hz,  Q’s  are  inferred  assuming  the  center  fre¬ 
quencies  of  0.9,  1.8  and  3.7  Hz. 

The  R~^'^  and  R~^'^  assumptions  bracket  the  average  crustal  Qp  of  270  for  this 
model.  There  is  no  real  significant  difference  between  the  earthquake  and  explosion 
determinations  of  Q  with  this  procedure. 


Pg  BEST  FIT  EXPONENTS 

EXPLOSION  EARTHQUAKE 


0.625-1.25  Hz 

-1.3 

-1.4 

1.25-2.5  Hz 

-2.5 

-2.6 

2.5-5.0  Hz 

-4.6 

-3.7 

APPARENT  Op„ 

EXPLOSION  EARTHQUAKEh=8km 


EXPONENT 

f=0.9 

1.8 

3.7 

f=0.9 

1.8 

3.7 

-1 

660 

317 

255 

507 

276 

338 

-5/6 

447 

280 

242 

378 

248 

317 

-1/2 

273 

233 

223 

245 

209 

280 

Contrary  to  the  suggestion  by  Haskell  (1966)  that  Pg  will  exhibit  a  leakage 
attenuation,  we  find  that  the  apparent  Q’s  for  the  synthetics  are  indistinguishable  from 
the  average  compressional  Q  of  the  crust.  Haskell  pointed  out  that  Pg  may  be  inter¬ 
preted  as  a  superposition  of  crustal  leaking  modes,  and  therefore  would  have  an 
apparent  spatial  attenuation  greater  than  the  intrinsic  attenuation  of  the  crust.  It  would 
seem  that  at  least  for  models  such  as  BR2,  that  the  leakage  attenuation  of  Pg  is  small. 
Haskell’s  analysis  was  motivated  in  part  by  the  Press’s  (1966)  observation’s  of  Qpg  = 
260  while  Qig  =  450.  Press’s  observations  for  Q^g  are  clearly  out  of  line  with  the 
majority  of  other  estimates  listed  in  Table  3.2,  although  his  Qpg  observations  are  con¬ 
sistent  with  those  of  Taylor,  et  al.  (1988).  It  is  unclear  why  Press’s  estimate  for  cru¬ 
stal  shear  Q  is  high  with  respect  to  the  other  estimates  listed  in  Table  3.2. 

The  same  approach  as  used  with  Pg  is  now  used  to  analyze  the  bandpass  filtered 
Pn  waveforms  and  leads  to  the  following  tabulated  results. 
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Pn  BEST  FIT  EXPONENTS 


EXPLOSION 

EARTHQUAKEh=8km 

0.625-1.25  Hz 

-2.9 

-2.8 

1.25-2.50  Hz 

-1.9 

-2.4 

2.50-5.0  Hz 

-3.4 

-3.2 

APPARENT  Qp, 

EXPLOSION _ EARTHQUAKE 


EXPONENT 

f=0.9 

1.8 

3.7Hz 

f=0.9 

1.8 

3.7Hz 

-1 

91 

385 

327 

97 

221 

338 

-5/6 

86 

327 

298 

88 

201 

317 

-1/2 

75 

252 

262 

77 

167 

275 

There  are  obviously  frequency  dependent  effects  in  the  Pn  Green’s  functions.  The  Q 
estimates  correspond  to  a  Q  proportional  to  and  for  exponents  -1, 

*5/6,  and  -1/2  respectively,  if  a  power  law  were  chosen  as  the  parameterization,  Q  = 
Qq^^.  The  effect  of  an  apparent  Qp^  nearly  proportional  to  frequency  is  profound  for 

—  dt*  R 

the  parameter  that  determines  the  shape  of  the  spectrum,  t*  =  t*  +  f  -rr  =  — ;r~ 

df  _  cQo 

(l-Ti)/"^.  For  the  spread  of  values  quoted  above  we  find  that  the  predicted  t*  at  500 
km  would  be  in  the  0.05  to  0.1  second  range.  This  is  within  the  wide  range  of  values 
de^mined  from  LRSM  stations  recording  HARDHAT  Pn  listed  in  Table  3-3.  If  such 
a  t*  were  interpreted  as  an  apparently  frequency  independent  Q,  it  would  imply  an 
apparent  Q  of  ^tween  650  and  1300  for  Pn.  The  estimates  are  also  consistent  with 
the  results  from  spectral  ratios  of  the  synthetics  presented  above.  The  values  in  the  1 
to  2  Hz  range  are  roughly  in  agreement  with  Taylor,  et  al.  ’s  (1988)  estimate  for  Pn 
Q’s  of  250  in  this  range.  The  best  fit  exponential  decays  for  aU  bands  are  in  the  R~^ 
to  R~^^  range  consistent  with  observations  of  recordings  with  similar  bandwidths. 

Finally,  the  analysis  for  Lg  was  performed  on  two  earthquake  sources  at  depths  of 
0.5  and  8  km  since  the  explosion  synthetics  had  weak  Lg  amplitudes  at  distances 
greater  than  4(X)  km.  For  the  1.25  to  2.5  and  2.5  to  5.0  Hz  bandwidths  the  Lg  ampli¬ 
tudes  at  large  distances  were  too  small  to  see  above  the  numerical  noise  and  Pg  coda. 
In  these  cases,  we  could  only  determine  upper  limits  to  the  exponents  and  to  the 
apparent  Q.  In  the  case  of  the  8  km  deep  earthquake,  the  Lg  signal  is  contaminated  by 
Sn  and  the  Q  is  higher  than  the  crustal  of  the  model.  Lg  was  difficult  to  measure  at 
high-firequencies  since  the  synthetic  signal  is  too  small  at  distances  greater  than  400 
km.  The  Lg  excitation  for  the  explosion  is  weak.  Lg  excitation  even  for  a  shallow 
earthquake  (h  =  0.5km)  is  weak  and  the  attenuation  imposes  a  severe  limitation  on  the 
observed  Lg  at  distance.  The  Lg  is  often  contaminated  by  Pg  coda.  Since  spatial 
attenuation  near  1  Hz  is  within  reasonable  agreement  with  observations  we  are  faced 
with  the  conclusion  that  crustal  Q  must  increase  with  frequency  in  the  basin  and  range 
or  we  would  not  be  able  to  see  Lg  at  distances  greater  than  300  km  at  fi^uencies 
above  1  Hz. 
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Le  BEST  FIT  EXPONENTS 

EARTHQUAKE(h=8km)  EARTHQUAKE(h=0.5km) 

0.625- 1.25Hz 

-2.8 

-3.5 

1.25  -2.5  Hz 

<-2.4 

- 

2.5  -  5.0  Hz 

<-3.2 

- 

APPARENT  Qr  „ 


EARTHQUAKE(h=8km) 

EARTHQUAKE(h=0.5km) 

EXPONENT 

f=0.9 

1.8 

3.7Hz 

f=0.9  Hz 

-1 

263 

<304 

<530 

156 

-5/6 

241 

<288 

<503 

145 

-1/2 

187 

<256 

<454 

129 

3.5  Conclusions 

Using  a  realistic  crust  and  upper  mantle  model  for  the  Basin  and  Range,  we  have 
demonstrated  that  spatial  attenuation  of  regional  phases  can  be  properly  modeled  near 
1  Hz  consistent  with  the  literature.  Synthetics  in  this  model  predict  no  significant 
differences  between  the  apparent  spatial  attenuation  rates  between  earthquakes  and 
explosions.  The  velocity  and  Q  model  will  be  used  in  a  subsequent  section  of  this 
report  to  propagate  seismic  waves  from  more  complicated  explosion  models  to  regional 
distances.  Although  our  main  goal  was  to  validate  a  model  to  be  used  to  simulate 
regional  propagation  suitable  for  discrimination  studies,  we  have  found  several  interest¬ 
ing  ancillary  results. 

First,  the  Airy  phase  cylindrical  spreading  approximation  at  1  Hz  should  yield  an 
excellent  measure  of  the  "average"  shear  Q  in  the  crust.  Sinoilar  results  were  found  for 
the  apparent  attenuation  and  geometrical  spre^ng  of  synthetic  Pg. 

Second,  the  frequency  dependence  of  synthetic  Pn  mimics  a  Q  that  increases  with 
increasing  frequency.  This  was  found  despite  the  use  of  a  frequency  independent 
attenuation  in  the  model.  The  functional  dependence  of  this  Pn  propagation  will  surely 
depend  on  the  fine  points  of  the  velocity  model  and  should  be  more  closely  investi¬ 
gated. 

Third,  the  explosion  point  source  is  a  very  poor  source  of  Lg.  This  has  been 
noted  before,  and  various  mechanisms  have  been  proposed  for  the  generation  of  SV, 
including  P  to  SV  conversion  at  interfaces  (Bennett,  et  al.  1987),  more  general 
scattering  (Gupta  and  Blandford,  1983),  and  spall  (Patton,  1988).  It  is  clear  that,  the 
problem  becomes  even  worse  when  realistic  attenuation  is  introduced  into  the  regional 
Green’s  function. 

For  both  the  earthquake  and  the  explosion  generated  synthetic  Lg,  the  high  fre¬ 
quency  amplitudes  are  too  low  to  match  observations  if  a  frequency  independent 
attenuation  operator  is  used.  Clearly,  this  suggests  that  the  empirical  frequency  depen¬ 
dence  of  the  Lg  Q  between  and  actually  reflect  the  frequency  dependence  of 
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the  average  shear  Q  of  the  crust.  Additional  study  is  required  to  determine  if  Q(f) 
models  with  more  complicated  depth  dependence  will  require  that  the  entire  crust 
reflect  such  a  dependence  or  whether  the  phenomenon  is  isolated  to  some  specific 
depth  range. 
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IV.  TWO-DIMENSIONAL  SOURCE  FAR-FIELD  P  AND  SV  WAVEFORMS 

4.1  Introduction 

In  this  section  we  examine  the  far-field  P-  and  SV-waves  radiated  by  several  non¬ 
linear  axisymmetric  finite  difference  simulations.  There  are  "Shagan  River"  simulations 
at  four  different  depths  of  burial  and  constant  yield.  Simulations  for  a  Pahute  Mesa 
tuff  model  were  conducted  at  three  different  depths  of  burial  and  constant  yield.  And 
finally,  there  are  three  simulations  intended  to  model  explosions  in  prestrcssed  granite. 

Although  our  principle  aim  is  to  compute  regional  synthetics,  there  are  several 
reasons  for  studying  the  far-field  body  waves  from  these  simulations.  The  body  waves 
are  simpler  and  can  be  compared  directly  to  observed  waveforms.  The  simpler 
Green’s  function  facilitates  a  more  direct  study  of  the  physics  of  the  source. 

The  Shagan  River  simulations  have  been  previously  discussed  by  Day,  et  al. 
(1986b)  in  order  to  derive  magnitude  corrections  for  the  January  15,  1965  presumed 
cratering  explosion  at  Shagan  River.  The  predicted  waveform  statistics  of  "a",  "b", 
and  "c"  phase  were  compared  to  observations  in  order  to  correct  the  presumed  crater¬ 
ing  WWSSN  mb  to  an  equivalent  contained  explosion.  We  re-examine  this  analysis 
and  estimate  the  150  kt  mb  magnitude  limit  at  this  test  site. 

Similarly,  we  examine  the  P  waveform  statistics  for  the  Pahute  Mesa  tuff  model. 
In  light  of  the  observations  of  Taylor,  et  al.  (1988)  that  normal  buried  events  at  NTS 
are  relatively  low-frequency  with  respect  to  overburied  events,  we  examine  the  P  and 
SV  waves  radiated  by  the  Pahute  Mesa  tuff  models  three  depths  of  burial. 

Finally,  the  body  waves  from  three  calculations  of  explosions  in  prestressed  gran¬ 
ite  are  examined.  These  "PILEDRIVER"  simulations  are  examined  to  observe  the 
effects  of  tectonic  release  (prestress)  on  the  short-period  P-SV  waveforms.  These 
simulations  were  previously  used  by  Day,  et  al.  (1986a)  to  demonstrate  the  effects  of 
tectonic  release  on  long-period  Rayleigh  waves.  We  extend  the  seismic  analysis  to  the 
effects  of  tectonic  release  on  far-field  broadband  P  and  SV  body  waves  from  explo¬ 
sions  in  granite. 

4.2  Shagan  River  Simulations 
4.2.1  Waveform  Statistics 

McLaughlin,  et  al.  (1985  and  1986)  presented  evidence  that  statistical  differences 
could  be  seen  in  teleseismic  P-waves  from  Shagan  River  explosions  diagnostic  of  the 
depth  of  burial.  The  statistics  of  the  relative  amplitudes  of  the  first  positive  swing, 
"a",  the  first  negative  swing,  "b",  and  the  second  positive  swing,  "c",  or  "max"  of  the 
WWSCN  P-wave  were  examined  (see  Figure  4.1).  It  was  found  that  the  network 
statistics  for  Log(Pmax/Pa),  Log(Pb/Pa),  and  Log(Pmax/Pb)  showed  that  the  January 
15,  1965  presumed  cratering  explosion  at  Shagan  River  had  an  magnitude  0.10  to 
0.30  units  less  than  contained  explosions  with  comparable  "a"  phase  amplitude.  The 
yield  of  the  January  15,  1965  event  is  presumed  to  be  125  kt  (Marshall,  et  al.,  1979; 
Rodean,  1979).  Therefore,  this  suggests  an  alternative  manner  by  which  the  Shagan 
River  test  site  could  be  calibrated.  Day,  et  al.  (1986b)  used  axisymmetric  nonlinear 
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d  =  Maximujn 


d 


Definition  of  Pa,  Pb,  and  Pmax  amplitudes  read  from 
film  records.  Pb  and  Pmax  are  1/2  the  peak-to-peak  b 
and  max  measurements  while  Pa  is  the  zero-to-peak 
amplitude  measurement. 


Figure  4.1. 
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finite  difference  simulations  to  predict  teleseismic  P-waveforms  for  Shagan  River 
events  of  125  kt  at  several  depths  of  burial.  TTieir  results  were  in  apparent  agreement 
with  McLaughlin,  et  al.  However,  as  we  will  see,  the  predicted  Log(Pmax/Pa)  from 
the  synthetics  is  not  in  agreement  with  observations  for  contained  shots. 

In  an  effort  to  understand  this  discrepancy,  we  have  re-examined  the  waveform 
statistics  for  the  S-CUBED  nonlinear  simulations.  Furthermore,  we  have  performed 
the  analysis  on  simulations  for  the  Pahute  Mesa  tuff  model.  The  analyses  have  been 
repeated  with  different  attenuation  models,  with  both  constant  t*  and  frequency  depen¬ 
dent  t*(f).  The  effects  of  typical  crustal  layering  at  the  receiver  have  also  been  exam¬ 
ined  for  the  Pahute  Mesa  synthetics. 

Figure  4.2  shows  predicted  teleseismic  P-wave  WWSSN  waveforms  for  the 
Shagan  River  simulations  of  125  kt,  at  depths  of  200  meters  (optimal  cratering  depth), 
680  meters  (normal  containment  depth)  and  980  meters  (overburied).  The  seismograms 
include  the  effects  of  attenuation  (constant  t*  =  0.7  sec)  and  instrument  response  for 
three  take  off  angles  of  10,  20  and  30  degrees.  We  are  concerned  here  with  the 
equivalent  teleseismic  signals  represented  by  the  synthetics  for  a  takeoff  angle  of  10 
degrees.  We  see  that  the  waveform  changes  morphology  as  depth  is  increased;  the  "c" 
phase  becomes  larger  relative  to  the  "a"  phase  as  depth  is  increased.  The  "a"  phase 
amplitude  is  largely  determined  by  the  amplitude  of  the  direct  P  wave,  the  attenuation 
operator  and  the  instrument  response.  The  later  "b"  and  "c"  amplitudes  are  affected  by 
the  constructive  and  destructive  interference  of  the  direct  and  later  portions  of  the  P- 
wave  as  well  as  the  attenuation  and  instrument  operators. 

Figure  4.3a  shows  the  far-field  displacement  pulses  without  the  attenuation  and 
instrument  operators  included  in  Figure  4.2.  From  these  waveforms  we  see  that  the 
simulations  all  produce  a  negative  polarity  "pP"  but  that  the  amplitude  and  timing  is  a 
function  of  depth.  The  deeper  explosions  show  negative  pulses  somewhat  attenuated 
with  respect  to  the  direct  P  wave  and  perhaps  distorted  by  propagation  through  the 
nonlinear  region  above  the  explosion.  The  negative  phase  for  the  cratering  simulation 
is  related  to  crater  formation  and  is  only  referred  to  as  a  "pP"  to  simplify  terminology. 
These  pulses  are  compared  in  Figure  4.3b  with  deconvolutions  from  Der,  et  al.  (1987) 
for  the  January  15,  1965  event  and  several  other  Shagan  River  explosions  of  similar 
magnitude.  We  see  that  the  cratering  simulation  is  in  fair  agreement  with  the  decon¬ 
volution  of  the  January  15,  1965  event  while  the  680  meter  simulation  does  not  have 
the  large  negative  polarity  pulse  prominent  in  the  deconvolutions  of  each  of  the  5.8 
or  5.9  events.  The  680  meter  simulation  appears  to  show  a  smaller  and  delayed  "pP" 
when  compared  to  the  deconvolutions. 

Figure  4.4  shows  the  predicted  statistics  for  Log(Pmax/Pa)  versus  Log(Pmax/Pb) 
from  the  simulations  shown  in  Figure  4.2  compared  to  some  observed  values 
(McLaughlin,  et  al. ,  1986).  The  observed  data  is  only  for  events  with  mi,  (Pmax)  > 
5.5.  The  simulations  predict  Log(Pb/Pa)  to  be  about  0.3  which  is  in  general  agreement 
with  observations  which  range  from  0.25  to  0.50.  The  observations  show  a  distribu¬ 
tion  of  Log(Pmax/Pa)  and  lx)g(Pmax/Pb)  with  the  presumed  cratering  shot  with 
minimum  Pmax  amplitude  relative  to  "a"  or  "b"  amplitude.  'I  he  simulations  predict 
Log(Pmax/Pa)  and  Log(Pmax/Pb)  for  the  "normal"  contained  shot  nearly  identical  with 


26 


18780611  5.9 


SSS-R-88-9844 


u  bC  C 
‘.2  C  aS 

5  2  ‘i: 

®  ^  S 

(O  u 

•t;  a 

?  2  e  - 

£  w 

>  ir 

^  M.  QJ  W 
“ 

bC  jN|  lO  J5 

Ih  ^  — 

fl;  ^ 

j3  O  4) 

OD  r.  lO  > 
4)  H  JO  Rj 

S;  ® 

*  -H 

^  »  s 

•  »«  ei 

0^  .2 

0)  r3 
’C  c  4> 
C  D. 

■  •S  *4^0* 

*7  o  -r  r3 

i  oo  ? 

aj  to  « 


29 


SHAGAN  RIVLR  Pmax.Pb.Rcj  Statistics  NTS  Pmax.Pb.Po  Statistics 


O  00  CD  CM 

d  o  c  d 


(od/xoujd)001 


O  CO  CD  CM 

d  d  d  d 


(0^/xDUJd)001 


30 


Figure  4.4.  Observed  WWSSN  waveform  statistics  Log(Pmax/Pa)  and  Log(Pmax/Pb)  compared 

to  synthetic  predictions.  The  synthetic  Log(Pb/Pa)  predictions  are  approximately 
correct  but  the  synthetics  under  predicted  Pmax  amplitudes  relative  to  the  Pa  and 
Pb  amplitudes. 
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the  observed  values  for  the  cratering  shot.  Furthermore,  they  predict  that 
Log(Pmax/Pb)  =  0.0  for  the  200  meter  depth  shot,  in  contrast  to  observations  for  the 
cratering  shot  of  about  0.15.  Clearly,  the  simulations  are  under  predicting  the  Pmax  in 
all  cases  although  they  do  an  acceptable  job  of  predicting  the  relative  amplitudes  of  Pa 
and  Pb. 

We  are  left  with  the  conclusion  that  the  nonlinear  contained  simulations  under 
predict  the  "pP"  amplitude  that  is  responsible  for  the  interference  effect  that  produces 
the  large  "c"  phase  amplitude.  This  is  based  on  the  relative  Log-amplitude  statistics  as 
well  as  the  deconvolution  results  that  generally  show  distinct  pP  arrivals  for  Shagan 
River  explosions  in  the  mi,  5.5  to  6.2  range. 

We  see  that  this  problem  also  exists  for  the  Pahute  Mesa  tuff  model.  Figure  4.4 
shows  the  Log(Pmax/Pa)  and  Log(Pmax/Pb)  predictions  for  NTS  based  on  the  Pahute 
Mesa  tuff  calculations  compared  to  observations  for  three  Pahute  Mesa  events. 
Again,  the  value  of  Log(Pb/Pa)  =  0.2  is  roughly  in  agreement  with  the  data  (about 
0.3),  but  the  Pmax  amplitude  is  under  predicted  by  about  a  factor  of  2  relative  to  the 
Pa  and  Pb  amplitudes.  We  will  come  back  to  this  in  a  later  section. 

4.2.2  Implications  for  Magnitude:YieId  at  Shagan  River 

Given  that  the  simulations  incorrectly  predict  the  Pmax  phase,  but  appear  to  give 
agreement  on  Pa  and  Pb,  we  re-examine  the  calibration  of  the  Shagan  River  test  site 
using  the  presumed  cratering  event  as  a  125  kt  reference  event.  From  the  simulations, 
we  have  that  Log(Pa)contained  -  Log(Pa)cratering  =  -0.128,  and  Log(Pb)contained  - 
Log(Pb)cratering  =  0.0.  As  pointed  out  by  Day,  et  al. ,  the  2D  synthetics  predict  an 
"a"  phase  coupling  scaling  with  depth  nearly  equal  to  that  predicted  by  Murphy 
(1977);  coupling  inversely  proportional  to  the  cube-root  of  depth.  Such  scaling  would 
predict  the  shallower  event  to  have  a  larger  RDP  by  Log(680/200)/3  =  0.177  units. 
Since  the  event  cratered  in  the  first  few  tenths  of  a  second,  the  simple  depth  scaling 
should  only  be  expected  to  give  a  qualitative  estimate  of  the  depth  scaling  for  the  ini¬ 
tial  phase  and  is  expected  to  progressively  fail  for  later  times  in  the  P-waveform. 
From  the  synthetics,  we  expect  that  a  contained  explosion  with  a  yield  of  125  kt  would 
have  had  an  mi,(Pa)  0.128  less  than  the  January  15,  1965  event  but  with  the  same 
mt(Pb).  Indeed,  if  we  examine  the  magnitudes  from  in  Table  4.1,  we  see  that  the 
events  of  November  30,  1969  and  September,  15,  1978  have  /^^(Pb)  values  close  to 
the  presumed  cratering  explosion  and  that  their  (Pb)-/n^  (Pa)  differences  (5.77- 
5.38=0.39,  and  5.67-5.35=0.32)  are  close  to  what  is  predicted  from  an  adjustment  of 
the  January  15,  1965  cratering  explosion  mi,  (Pa)-mi,  (Pa),  5.72-5.48+0.128=0.37.  Con¬ 
sequently,  we  estimate  that  for  the  WWSSN  magnitude  estimation  procedure  used  by 
McLaughlin,  et  al.  (1986)  that  a  normally  contained  125  kt  explosion  corresponds  to 
m^(Pa)  =  5.48(0.05)  -  0.13  =  5.35(0.05)  and  m^(Pb)  =  5.72(0.05).  In  Appendix  A  we 
use  these  corrected  magnitudes  to  estimate  the  150  kt  limit  at  the  Shagan  River  test 
site. 
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TABLE  4.1  MLE  MAGNITUDES  (McLaughlin  et  al.  (1986)) 
YY/MM/DD _ m/j  (Pa) _ m/,  (Pb) _ nii,  (Pmax) 


65/01/15 

5.48(0.05) 

5.72(0.05) 

5.87(0.04) 

68/06/19 

4.59(0.05) 

5.02(0.04) 

5.27(0.04) 

69/11/30 

5.38(0.04) 

5.77(0.03) 

5.97(0.03) 

73/07/23 

5.76(0.05) 

6.02(0.05) 

6.20(0.04) 

73/12/14 

5.16(0.05) 

5.49(0.05) 

5.74(0.05) 

75/04/27 

4.99(0.12) 

5.32(0.10) 

5.53(0.10) 

76/07/04 

5.24(0.08) 

5.61(0.05) 

4.91(0.06) 

76/12/07 

4.93(0.15) 

5.42(0.12) 

5.61(0.13) 

78/06/11 

5.25(0.05) 

5.57(0.05) 

5.85(0.05) 

78/09/15 

5.35(0.06) 

5.67(0.05) 

5.86(0.06) 

79/06/23 

5.72(0.06) 

5.96(0.06) 

6.18(0.06) 

80/09/14 

5.58(0.09) 

5.90(0.08) 

6.14(0.08) 

4.3  Pahute  Mesa  Simulations 

In  the  following  section,  we  describe  seismic  waveforms  from  simulations  done  in 
a  layered  medium  which  models  explosion  emplacement  at  Pahute  Mesa.  The  model 
has  three  layers  representing  a  very  low  velocity  surface  layer,  and  a  strong  velocity 
increase  at  the  water  table.  As  with  the  Shagan  River  calculations,  the  simulations  for 
Pahute  Mesa  were  done  for  a  yield  of  125  kt,  and  at  a  cratering  depth  (200  m),  an 
optimal  depth  (680  m),  and  an  overburied  depth  (980  m).  The  2D  calculations  include 
the  full  nonlinear  response  of  the  material  as  well  as  the  effects  of  gravity,  cracking, 
spall,  and  the  ballistic  material  response  in  the  case  of  the  cratering  calculation  (200 
m).  A  "generic"  Pahute  Mesa  structure  was  constructed  based  on  properties  from  a 
number  of  drill  holes  in  the  area.  The  seismic  velocity  model  at  the  source  is  given  in 
Table  4.2.  Other  material  properties  of  the  model  are  given  by  Day,  et  al.  (1986). 


TABLE  4.2  PAHUTE  MESA  TUFF  MODEL 


layer  thickness  (km) 

P  (km/s) 

a  (km/s) 

p  (gm/cc) 

depth  (km) 

0.122 

0.661 

1.204 

1.00 

0.0-0.122 

0.460 

1.109 

2.025 

1.95 

0.122-0.582 

OO 

1.581 

2.887 

2.00 

0.582-OO 

The  contained  calculations  were  constrained  so  as  to  match  the  cavity  sizes  for 
events  of  these  depths  and  yield.  It  rs  interesting  to  note  that  the  calculations  for  the 
125  kt,  980  meter  depth  simulation  were  subsequently  found  to  match  the  surface 
ground-zero  peak  vertical  velocity  and  the  free-fall  (spall)  observed  time  for  SCOTCH 
to  within  10%.  Tensile  failure  occurred  in  the  simulations  at  a  number  of  horizons 
above  the  working  point,  but  the  actual  spall  in  the  calculations,  defined  as  a  period  of 
free-fall  (-1  g),  was  found  to  be  confined  to  the  weak  surface  layer  (upper  122  m). 
Events  of  this  size  and  depth  at  Pahute  Mesa,  for  which  data  is  available,  show  spall 
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reaching  to  depths  as  much  as  1/3  the  depth  of  burial  (H.  Patton,  personal  communica¬ 
tion).  Detailed  ground  motion  records  above  the  working  point  could  serve  as  impor¬ 
tant  constraints  for  future  modeling  efforts,  particularly  if  spall  openings  and  closures 
are  important  to  regional  phase  excitation. 

4.3.1  Waveform  Statistics 

We  have  computed  the  far-field  signal  along  rays  leaving  the  bottom  of  the 
source  structure  for  three  take-off  angles  (10,  20  and  30  degrees).  These  three  angles 
correspond  to  a  teleseismic  P  wave  at  60  degrees  (17  km/s),  Pn  apparent  velocity  in 
the  9  to  10  degree  range  (8.5  km/s),  and  a  crustal  Pg  apparent  velocity  (5.8  km/s).  In 
addition,  for  comparison  with  the  waveforms  of  observed  seismograms,  the  signals 
were  convolved  with  various  attenuation  operators  and  a  WWSSN  instrument  response. 

The  dependence  of  the  unfiltered  signals  on  take-off  angle  and  depth  of  burial  is 
shown  in  Figures  4.5a  and  b,  respectively.  The  first  peak,  the  direct  arrival  P,  is 
unaffected  within  this  range  of  take-off  angles,  but  the  trailing  signal,  the  free  surface 
interaction  phase  pP,  is  affected.  At  each  depth,  pP  is  largest  and  is  delayed  most  at 
the  steepest  take-off  angle  (10®  ).  This  is  predicted  by  elastic  theory.  The  increase  in 
pP  amplitude  is  greatest  at  the  deepest  depth.  The  direct  arrivals  at  depths  of  200  and 
680  m  are  nearly  the  same  amplitude,  but  the  amplitude  from  a  depth  of  980  m  is 
diminished  to  about  40%  of  the  others.  This  is  due  to  differences  in  material  proper¬ 
ties  with  depth.  The  amplitude  of  pP  relative  to  P  is  largest  for  the  deeper  shots, 
showing  the  effect  of  containing  the  explosion. 

To  isolate  the  effects  of  nonlinear  interactions  with  the  free  surface,  we  have 
compared  the  2D  simulations  discussed  above  with  comparable  ID  source  simulations. 
The  ID  solutions  were  obtained  by  first  performing  a  nonlinear  calculation  in  a  whole 
space  whose  material  properties  were  the  same  as  those  at  the  depth  of  burial  in  the 
2D  calculation,  and  then  extracting  the  reduced  velocity  potential  (RVP).  The  RVP 
was  convolved  with  the  elastic  source  structure  transfer  function  to  obtain  the  far-field 
wave  leaving  the  source  structure,  which  can  be  compared  directly  to  the  comparable 
2D  results.  The  2D  and  ID  waveforms  are  overlain  in  Figure  4.6. 

As  might  be  expected,  the  difference  between  the  ID  and  2D  free  surface  interac¬ 
tion  is  most  apparent  at  the  shallowest  depth  (200m),  where  cratering  occurred.  The 
direct  waves  are  the  same,  but  the  large  negative  peak  predicted  by  an  elastic  free  sur¬ 
face  reflection  is  almost  absent  in  the  2D  simulation.  This  is  the  case  for  all  three 
take-off  angles.  At  deeper  depths,  the  nonlinear  pP  is  generally  delayed  by  0.2 
seconds  relative  to  the  elastic  pP.  The  reduction  in  amplitude  of  pP  due  to  the  non¬ 
linear  free  surface  interaction  is  greatest  (  ~  1/3  )  at  the  shallowest  take-off  angle  ( 
30®  ),  and  least  (  25  %  )  at  the  steepest  take-off  angle  (  10®  ). 

To  compare  these  waveforms  with  features  of  observed  seismograms,  we  con¬ 
volved  the  time  series  discussed  above  with  an  attenuation  operator  and  an  instrument 
response.  Although  we  examine  the  effects  of  other  attenuation  operators  later  in  this 
paper,  we  begin  with  a  frequency-independent  operator  with  t*  =0.5.  A  short-period 
WWSSN  instrument  response  was  used.  The  synthetic  seismograms  for  the  three 
take-off  angles  are  overlain  at  the  three  depth  of  burials  in  Figure  4.7.  The  first  two 
peaks,  the  "a"  and  "b"  phases  are  the  same  for  each  take-off  angle  at  all  depth  of 
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burials,  while  the  "c"  phase  changes  slightly  in  amplitude  and  period  as  a  function  of 
take-off  angle  for  the  cratering  shot.  The  increase  in  the  ratio  of  pP  to  P  is  most 
apparent  at  the  siecpest  take-off  angle  (  10®  ).  Note  that  the  effect  on  die  ratio  of  "c" 
to  "b"  amplitude  is  small  at  all  take-off  angles.  The  effects  of  a  WWSSN  instrument 
response  and  t*  =  0.5  sec  are  shown  in  Figure  4.8. 

For  the  seismograms  discussed  above,  a  frequency-independent  Q  model  with 
t*  =0.5  was  used.  The  effects  of  using  the  frequency  dependent  Q  models  of  Der  and 
Lees  (1985)  can  be  seen  in  Figures  4.9.  The  models  are  designated  TT  for  a  tectonic- 
to-tectonic  path,  TS  for  a  tectonic-to- shield  path,  and  SS  for  a  shield-to-shield  path. 
Model  TT  is  most  attenuative,  SS  is  least  attenuative  and  TS  is  intermediate.  The  t*  at 
1  Hz  for  the  three  models  are  0.35  0.55  and  0.75  sec  for  SS,  TS,  and  TT  respectively. 
We  expect  that  the  paths  from  NTS  to  the  WWSSN  stations  are  a  mix  of  paths  brack¬ 
eting  the  TS  and  TT  models. 

In  Figure  4.9,  the  waveform  plots  are  scaled  by  the  largest  of  the  traces,  so  that 
the  relative  amplitudes  are  apparent.  As  expected,  the  SS  model  is  largest.  The 
seismograms  for  the  TS  model  and  the  frequency-independent  model  track  each  other 
closely.  The  differences  between  the  SS  model  and  the  others  is  apparent  in  the 
seismograms  for  the  cratering  shot,  where  the  inflection  in  the  "c"  phase,  due  to 
separation  of  P  and  pP,  can  be  seen  for  SS  but  shows  up  as  a  widened  peak  for  the 
other  Q  models.  The  separation  of  P  and  pP  is  apparent  to  different  degrees  for  the 
four  Q  models  in  the  optimally  buried  (680  m)  shot. 

As  was  done  in  Day,  et  al.  (1986b),  measurements  of  the  "a",  "b"  and  "c"  phases 
and  their  relative  values  were  made  for  the  synthetic  seismograms.  Their  values  are 
listed  in  Tables  B.la  through  B.le  of  Appendix  B.  Several  additional  figures  are 

shown  in  Appendix  B.  The  ratios  of  the  "b"  and  "c"  phases  (  peak-to-peak)  to  the 

"a"  (zero-to-peak)  are  shown.  These  measurements  are  referred  to  as  ^  bc^  and 

P  ah  P he 

Pa,  respectively.  These  ratios,  and  generally  less  than  0.25  (a  factor 

of  1.8)  for  the  various  models  for  both  the  contained  and  optimally  buried  shots,  with 
the  contained  somewhat  larger.  Observations  are  typically  about  0.6  (a  factor  of  4), 
indicating  that  the  simulated  pP  is  delayed  and/or  is  too  small  compared  to  observa¬ 
tions.  As  can  be  seen  from  Figure  4.9,  and  examination  of  Tables  B.la  through  B.le, 
the  Log(Pa/Pb)  and  Log(Pc/Pa)  are  little  affected  by  the  choice  of  attenuation  model. 
The  scatter  due  to  the  most  divergent  attenuation  models  is  less  than  0.05  log  units. 

4.3.2  Far>Field  P-Wave  Spectra  and  Depth  Scaling 

Figure  4.10  shows  the  P-wave  power  spectra  for  the  three  simulations  appropriate 
for  a  10  degree  takeoff  angle.  Note  that  the  200  meter,  cratering  simulation,  has  an 
over  all  lower  level  above  0.5  Hz.  Because  we  are  interested  in  discrimination  prob¬ 
lems  associated  with  over  burial,  we  will  concentrate  on  the  two  680  and  980  meter 
depth  of  burial  simulations.  The  680  and  980  meter  simulations  are  rather  similar 
although  the  980  meter  spectrum  has  deeper  spectral  nulls.  In  order  to  bring  out  the 
differences  between  the  two  sets  of  spectra,  we  computed  spectral  ratios  for  the  three 
different  take-off  angles  averaged  over  several  bandwidths  (Figure  4.11).  We  see  that 
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2D  synthetic  P-waves  at  takeoff  angle  10*,  with  WWSSN  instrument,  and  four 
different  attenuation  operators.  All  traces  shown  at  same  scale. 


SYNTHETIC  TELESEISMIC  P-WAVE  SPECTRA 

PAHUTE  MESA  TUFF  MODEL,  125  Kl,  TAKEOFF  ANGLE  10 
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Figure  4.10.  Synthetic  far-field  P-wave  power  spectra  for  the  Pahute  Mesa  ti 

takeoff  angle  10*  and  depths  of  burial,  200,  680,  and  980  meters. 
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all  three  takeoff  angles  show  a  net  enhancement  of  the  680  meter  spectra  below  1  Hz, 
and  a  net  enhancement  of  the  980  meter  spectra  above  1  Hz.  The  net  difference  in 
Log(680/980)  below  i  Hz  and  above  2  Hz  is  about  0.15  log  units  or  a  factor  of  40%. 
This  low/high  frequency  difference  for  the  two  scaled  depths  is  only  visible  with  the 
spectral  averaging.  The  difference  is  in  agreement  with  what  is  predicted  by  Murphy 
(1977)  depth  scaling  in  tuff.  Figure  4.11  shows  the  predicted  spectral  ratio  for  these 
two  scaled  depths  with  the  appropriate  material  properties.  The  simple  depth  scaling 
theory  predicts  the  low/high  frequency  difference  properly  but  not  the  actual  level  of 
the  ratio. 

Taylor,  et  al.  (1988)  have  observed  that  overburied  explosions  at  NTS  are  rela¬ 
tively  enriched  in  high  frequencies  with  respect  to  adjacent  normal  buried  explosions. 
In  fact  they  have  found  that  these  overburied  events  fail  the  high/low  frequency 
discriminant  and  appear  to  look  like  earthquakes.  However,  the  small  explosions  in 
their  data  base  (m^  <  4)  that  fail  the  high/low  frequency  discriminant  are  100%  or 
more  overburied.  It  is  clear  however,  that  the  nonlinear  simulations  and  the  simple 
depth  scaling  theory  are  in  qualitative  agreement  with  the  observations. 

To  illustrate  the  consequences  of  the  depth  scaling,  we  use  Murphy’s  (1977) 
depth  scaling  theory  with  a  velocity  model  of  Leonard  and  Johnson  (1987).  The 
predicted  spectral  ratio  between  two  5  kt  explosions  at  depths  of  burial  200  and  500 
meters  are  shown  in  Figure  4.12.  A  factor  of  3  enhancement  is  predicted  for  the  high 
frequency  spectrum  of  the  overburied  explosion.  Roughly  half  of  this  shift  in  spec¬ 
trum  is  produced  by  the  shift  in  comer  frequency  due  to  the  increase  in  over  burden 
pressure  with  depth,  and  the  remaining  difference  is  due  to  the  differences  in  material 
.  properties  associated  with  the  increase  in  compressional  velocity  with  depth.  The  rela¬ 
tive  high  frequency  enhancement  with  respect  to  normal  depth  of  burial  increases  with 
decreasing  yield  and  with  increasing  over  burial  such  that  factors  of  10  are  predicted 
for  1  kt  explosions  overburied  100%  or  more  in  tuff  with  the  Leonard  and  Johnson 
velocity  model. 

4.3J  Far-Field  P-  and  SV-Wave  Spectra 

One  important  way  in  which  two-dimensional  sources  may  differ  from  one¬ 
dimensional  point  sources  is  the  radiation  of  direct  SV  energy.  A  complete  set  of  P 
and  SV  far-field  displacement  waveforms  are  shown  in  Figtues  4.13a,  b,  and  c  for 
take-off  angles  of  10,  20  ,  and  30  degrees  respectively.  In  each  case,  the  2D  simula¬ 
tions  for  200,  680,  and  980  meter  depth  are  shown  along  with  the  elastic  response  of 
the  model  convolved  with  a  ID  RDP  that  was  computed  for  the  680  meter  depth. 
Note  that  there  is  a  small  SV  arrival  on  the  nonlinear  2D  synthetics  that  does  not 
appear  in  the  ID  model.  This  is  the  direct  SV  produced  by  the  axisymmetric  explo¬ 
sion  model.  Upon  visual  inspection,  one  can  see  that  the  2D  SV  pulses  appear  broader 
(lower  frequency)  than  the  ID  elastic  model.  Also  the  peak  amplitude  of  the  2D  680 
meter  SV  synthetic  is  delayed  with  respect  to  the  elastic  propagator  model.  However, 
the  general  features  of  the  2D  and  ID  model  are  similar  and  the  peaks  and  troughs  of 
the  2D  SV  pulses  can  be  correlated  with  similar  peaks  and  troughs  of  the  ID  elastic 
model.  Consequently,  we  can  qualitatively  describe  the  differences  between  the  2D 
nonlinear  synthetics  and  the  ID  linear  model  as  if  the  elastic  P-SV  conversions  above 
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Figure  4.12.  Prediction  for  P-wave  spectral  ratio  between  two  hypothetical  5  kt  sources  at  200 
and  500  meter  depths  of  burial.  A  factor  of  3  high/low  frequency  difference  is 
predicted  based  on  cube  root  depth  scaling  and  a  typical  NTS  velocity /density 
model. 
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Figure  4.13a. 


P  (left)  and  SV  (right)  displacements  at  a  takeoff  angle 
of  10*  for  the  2D  sources  at  200  (top),  680,  and  980 
meters  and  the  ID  elastic  (bottom)  source.  All  traces 
shown  to  same  scale.  Note  the  direct  SV  in  the  2D 
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Figure  4.13b. 


P  (left)  and  SV  (right)  displacements  at  a  takeoff  angle 
of  20*  for  the  2D  sources  at  200  (top),  680,  and  980 
meters  and  the  ID  elastic  (bottom)  source.  All  traces 
shown  to  same  scale.  Note  the  direct  SV  in  the  2D 


sources. 
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Figure  4.13c.  P  (left)  and  SV  (right)  displacements  at  a  takeoff  angle 
of  30*  for  the  2D  sources  at  200  (top),  680,  and  980 
meters  and  the  ID  elastic  (bottom)  source.  All  traces 
shown  to  same  scale.  Note  the  direct  SV  in  the  2D 


sources. 
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the  working  point  were  attenuated  and  delayed.  The  2D  980  meter  depth  synthetic  is 
the  most  like  the  elastic  680  meter  depth  synthetic  in  appearance.  Judging  from  the 
timing  of  the  secondary  P  and  SV  pulses  it  appears  that  the  major  SV  source  from  the 
2D  nonlinear  axisymmetric  simulations  remains  the  conversion  of  P  to  SV  at  the  free 
surface  (pS).  The  attenuation  of  the  high  frequency  SV  above  the  source  increases 
with  increasing  take-off  angle  and  is  greatest  for  the  30  degree  take-off  angle. 

Figures  4.14a,  b,  and  c.  show  the  linear-linear  spectra  corresponding  to  the  time 
functions  of  Figures  4.13a,  b,  and  c  respectively.  Note  that  the  2D  nonlinear  sources 
have  greater  low  frequency  amplitude  and  less  high  frequency  amplitude  than  the 
elastic  simulation.  Also,  note  that  the  spectral  nulls  clearly  visible  in  the  ID  elastic  P 
and  SV  spectra  are  not  as  apparent  in  the  2D  nonlinear  spectra.  Clearly,  the  qualita¬ 
tive  explanation  that  pP  and  pS  are  attenuated  and  delayed  is  inadequate  to  explain  the 
spectral  character  of  the  nonlinear  simulations.  If  the  pP  and  pS  were  simply  linearly 
attenuated,  the  spectral  nulls  would  remain.  Furthermore,  the  low  frequency  ampli¬ 
tudes  have  been  increased  as  well  as  the  high  frequency  amplitudes  reduced.  The  non¬ 
linear  interactions  with  the  free  surface  and  the  interfaces  above  the  680  and  980  meter 
sources  are  responsible  for  phases  that  appear  in  the  time  domain  like  pP  and  pS  but 
have  spectral  character  unlike  the  direct  explosion  source. 

4.3.4  Broadband  Far<Field  P  Waveforms 

Figures  4.15a  and  b  compare  broadband  teleseismic  synthetics  for  the  Pahute 
Mesa  simulations  with  deconvolutions  by  Lyman,  et  al.  (1986)  and  Der,  et  al. 
(1987).  Recordings  from  the  EKA  array  were  deconvolved  of  instrument  and  an 
attenuation  operator  using  either  the  array  sum  (Lyman,  et  al.)  or  a  multichannel 
deconvolution  estimator  (Der,  et  al. ).  Lyman,  et  al.  removed  an  attenuation  operator 
with  t*  =  0.35  sec,  while  Der,  et  al.  used  a  t*  =  0.45  sec.  By  comparing  the  two 
deconvolutions  of  PURSE,  one  can  see  that  the  major  characteristics  of  the  waveform 
are  robust  with  respect  to  variations  in  attenuation  operator  and  deconvolution 
method. 

The  events  chosen  for  comparison  are  all  announced  by  DOE  to  be  in  the  20  to 
200  or  20-150  kt  ranges.  Their  maximum-likelihood  mj,’s  (Lyman,  et  al.,  1986)  are 
greater  than  5.5  and  depths  of  burial  between  599  meters  (PURSE)  and  978  meters 
(SCOTCH).  According  to  Springer  and  Kinneman  (1971),  SCOTCH  has  a  yield  of  150 
kt  and  therefore  a  scaled  depth  of  burial  of  180  meters/kr^^^. 

The  broadband  seismograms  (from  Lyman,  et  al.)  are  labeled  with  "ApP"  indi¬ 
cating  the  inferred  pP  arrival,  and  "As"  a  sometimes  significant  positive  arrival  follow¬ 
ing  the  "ApP".  In  the  case  of  PURSE  there  is  a  phase  labeled  "Aq"  which  Lyman,  et 
al.  suggest  may  be  aftershock  activity.  Hamilton,  et  al.  (1972)  report  that  PURSE 
reactivated  the  aftershock  zone  of  BENHAM  that  was  detonated  4.5  months  earlier. 
Lay  (1987)  and  Lay  and  Welc  (1987)  observe  no  more  significant  coda  complexity  for 
PURSE  than  for  other  Pahute  Mesa  explosions  as  measured  across  the  WWSSN  net¬ 
work.  Lay  (1987)  however,  claims  that  SCOTCH  exhibited  anomalous  P-coda  com¬ 
plexity  on  average  across  the  WWSSN  and  yet  shows  little  unusual  coda  in  Figure 
4.15b.  The  Der,  et  al.  deconvolution  for  the  somewhat  larger  explosion  (mle-m^, 
6.24),  FONTINA,  shows  coda  levels  similar  to  PURSE.  FONTINA  and  PURSE  were 
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Figure  4.14a.  P  (left)  and  SV  (right)  linear-linear  spectra  at  a  takeoff 
angle  of  10*  for  the  2D  sources  at  200  (top),  680,  and 
980  meters  and  the  ID  elastic  (bottom)  source.  All 
spectra  shown  to  same  scale.  Note  the  spectral  scalloping 
in  the  ID  source  and  the  depletion  of  high  frequencies  in 
the  2D  SV  spectra. 
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Figure  4.14b.  P  (left)  and  SV  (right)  linear-linear  spectra  at  a  takeoff 
angle  of  20*  for  the  2D  sources  at  200  (top),  680,  and 
980  meters  and  the  ID  elastic  (bottom)  source.  All 
spectra  shown  to  same  scale.  Note  the  spectral  scalloping 
in  the  ID  source  and  the  depletion  of  high  frequencies  in 
the  2D  SV  spectra. 
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Figure  4.14c.  P  (left)  and  SV  (right)  linear-linear  spectra  at  a  takeoff 
angle  of  30*  for  the  2D  sources  at  200  (top),  680,  and 
980  meters  and  the  ID  elastic  (bottom)  source.  All 
spectra  shown  to  same  scale.  Note  the  spectral  scalloping 
in  the  ID  source  and  the  depletion  of  high  frequencies  in 
the  2D  SV  spectra. 
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FIGURE  41.  P  SEISMOGRAMS  RECORDED  AT  ESKDALEWIR.  SCOTLAND  FROM  NTS 
EXPLOSION  BACKBEACH 


FI.PyRE  SEISMOGRAMS  RECORDED  AT  ESKDALEMUIR,  SCOTIAyP  FROM  NTS 
EXPLOSION  PANIR '  ~ 


Figure  4.15a.  Deconvolved  waveforms  from  Lyman,  tt  al.  for 
BACKBEACH  (top)  and  PANIR  (bottom). 
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Figure  4.15b.  Deconvolved  waveforms  from  Lyman,  et  al.  for  SCOTCH 
(top)  PURSE  (bottom),  and  from  Der,  et  al.  for  PURSE 
(middle).  The  two  deconvolved  waveforms  for  PURSE 
used  different  procedures  and  different  t*’s. 


53 


SSS-R-88-9844 


located  within  2  km  of  each  other  and  may  have  illuminated  some  nearby  scatterer  or 
FONTINA  may  also  have  re-activated  a  cluster  of  activity  that  was  observed  by  Ham¬ 
ilton,  et  al.  4  km  east  of  PURSE  and  would  have  been  within  2  km  of  FONTINA. 

The  980  meter  deep  2D  synthetic  waveform  (Figure  4.15c)  exhibits  a  1.1  sec 
pP-P  and  the  largest  positive  swing  in  the  "coda"  is  at  1.3  seconds  after  the  P  and  less 
than  1/2  the  amplitude  seen  in  the  deconvolutions.  The  apparent  pP  amplitude  of  the 
overburied  synthetic  is  substantially  less  than  most  of  the  observations. 

With  these  small  "pP"  signals  and  without  the  large  "As"  phases  in  the  synthetics 
it  is  not  possible  to  properly  simulate  the  Pmax/Pa  amplitude  ratios  observed  on  telese- 
ismic  WWSSN  seismograms.  These  secondary  arrivals  in  the  very  early  coda  of 
teleseismic  P  have  been  noted  by  Frasier  (1972)  and  many  others.  Given  that  this 
"phase"  is  so  prevalent  on  explosions  from  NTS  it  is  unlikely  that  it  is  a  scattered 
phase  but  rather  related  to  nonlinear  interaction  of  the  explosion  with  the  free  surface. 
Springer  (1974)  has  suggested  spall  closure  is  a  significant  mechanism  for  the  genera¬ 
tion  of  this  phase.  Although  the  2D  axisymtmtric  nonlinear  calculations  include  spall 
and  nonlinear  behavior  above  the  explosion,  they  may  underestimate  the  contribution 
of  spall  if  it  is  not  uniform. 

For  comparison  with  the  synthetics  at  shallower  take  off  angles,  Figure  4.16a 
shows  deconvolved  Pn  signals  from  a  set  of  events  at  the  broadband  stations  KNB  and 
MNV.  A  comparison  of  the  observations  with  synthetics  is  made  in  Figure  4.16b. 
The  agreement  between  synthetics  and  the  observations  is  acceptable  given  the 
observed  variation  between  the  two  stations  MNV  and  KNB.  The  most  striking  aspect 
of  both  the  observations  and  the  2D  synthetics,  is  the  lack  of  any  noticeable  "pP" 
phase.  This  is  in  contrast  to  the  elastic  Green’s  function  with  a  point  explosion  source 
which  predicts  a  prominent  pP  phase  (dashed  line  in  Figure  4.16b). 


TABLE  4.3  PAHUTE  MESA  EVENTS  USED  FOR  COMPARISON 
-ORDERED  BY  DEPTH  OF  BURIAL  (DOB)- 


EVENT 

nih 

Ms 

DOB(m) 

YIELD(kt) 

SCALED-DEPTH 

PURSE 

5.66 

4.74 

599 

20-200 

>  102 

BACKBEACH 

5.54 

- 

611 

20-150 

>  115 

FONDUTTA 

5.49 

- 

633 

20-150 

>  119 

PANIR 

5.66 

- 

681 

20-150 

>  128 

FARM 

5.63 

- 

689 

20-150 

>  130 

SCOTCH 

5.59 

4.45 

978 

150 

180 

4.4  Effect  of  Prestress  on  Far-Field  P  and  SV  Waves 

Day,  et  al.  (1986)  presented  nonlinear  models  of  tectonic  release  from  under¬ 
ground  explosions.  They  illustrated  their  analysis  with  synthetic  long  period  Rayleigh 
waves  from  a  suite  of  granite  source  calculations.  One  of  these  three  2D  calculations 
modeled  PILED RI\^R  (61  kt,  463  meter  depth  of  burial)  in  granite.  The  remaining 
two  calculations  were  performed  with  horizontal  prestress  (7.5  MPa  and  15Mpa)  to 
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Figure  4.15c.  Far-field  teleseismic  P-waves  (takeoff  angle  of  10*)  from 
the  2D  sources  at  200,  fi80  and  980  meter  depths  of 
burial  and  the  ID  RVP  at  680  meters.  Same  scale  as 
Figures  4.15a  and  b. 
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Figure  4.16a.  Deconvolved  Pn  waveforms  from  several  Pahute  Mesa 
events  recorded  at  KNB  and  MNV.  A  “pPn”  arrival  is 
not  generally  apparent  on  the  displacement  waveforms. 
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examine  the  effects  of  tectonic  release.  In  this  section,  we  examine  the  short-period 
far-field  body  waves  from  these  simulations  and  compare  the  no  prestress  and  prestress 
explosion  sources  with  a  ID  point  source  model  for  the  explosion. 

Figure  4.17a  shows  the  far-field  displacement  waveforms  for  the  three  2D  non¬ 
linear  calculations  and  a  ID  source  convolved  with  the  elastic  Green’s  function  for  a 
point  source  with  20°  takeoff  angle.  We  see  that  the  three  2D  calculations  are  very 
similar  with  an  apparent  "pP"  about  50%  the  amplitude  of  the  initial  P.  The  ID 
source  apparent  "pP"  is  larger  than  the  initial  positive  P  by  about  10%,  because  the 
ID  RDP  has  considerable  overshoot.  These  waveforms  are  convolved  with  an  instru¬ 
ment  response  and  an  attenuation  operator  (t*  =  0.8  sec)  in  Figure  4.17b.  We  see 
again  that  the  three  2D  P  waveforms  are  very  similar.  All  sources  give  the  same  Pa 
amplitude  and  that  the  ID  waveform  has  a  distinctively  larger  Pc  swing.  As  with  the 
Pahute  Mesa  and  Shagan  River  synthetic  s,  the  relative  amplitude  of  the  Pc  phase  to  the 
Pb  and  Pa  phase  is  smaller  than  the  ’  short-period  waveform.  In  order  to 
increase  the  Pc  phase  amplitude,  the  2D  sii*  Mlit.ons  would  be  required  to  have  larger 
"pP"  arrivals. 

The  far-field  SV  waveforms  for  three  takeoff  angles  (10°,  20°,  and  30°)  are 

shown  in  Figures  4.18a,  b,  and  c.  The  ID  source  has  the  largest  broadband  peak  SV 

amplitude  at  all  takeoff  angles.  The  2D  simulations  generate  a  small  direct  SV  wave 
that  precedes  the  P  to  S  conversion  in  the  elastic  Green’s  function.  The  differences 
between  the  2D  sources  are  smallest  for  10°  takeoff  angle  and  greatest  for  30°  takeoff 
angle.  The  direct  SV  is  largest  for  the  largest  prestress. 

The  broadband  displacement  waveforms  are  summarized  (takeoff  angle  30°)  in 
Figure  4.19a.  They  are  shown  all  at  the  same  scale  and  we  can  see  that  the  2D  SV 

pulses  are  generally  broader  than  the  2D  P  pulses  unlike  the  ID  SV  and  P  pulses 

which  have  similar  pulse  widths.  The  spectra  of  these  broadband  signals  are  shown  in 
Figure  4.19b.  We  see  that  the  2D  P  spectra  are  broader  than  the  ID  P  spectrum,  and 
that  the  2D  SV  spectra  are  more  sharply  peaked  than  the  ID  SV  spectrum.  As  is  evi¬ 
dent  in  the  time  domain  pulses,  the  2D  SV  spectra  have  a  dominant  peak  at  lower  fre¬ 
quency  than  the  2D  P  spectra. 

We  can  see  from  the  2D  body  waves  that  the  horizontal  prestress  has  little  effect 
on  the  amplitudes  and  waveforms  of  the  far-field  short  period  P  and  SV  signals.  The 
2D  P  waveforms  are  quite  different  than  the  ID  waveforms,  but  as  with  the  Pahute 
Mesa  and  Shagan  River  simulations,  the  Pc  phase  relative  to  the  Pa  and  Pb  phases 
appears  to  be  smaller  than  observed. 

4.6  Conclu.sions 

The  discrepancies  between  the  Pahute  Mesa  and  Shagan  River  synthetic  and 
observed  teleseismic  waveform  statistics,  together  with  the  deconvolution  results  indi¬ 
cate  that  the  nonlinear  free-surface  interaction  phase,  "pP"  is  too  small  in  the  synthet¬ 
ics.  It  is  possible  that  this  is  related  to  the  fact  that  spall  was  observed  in  the  synthet¬ 
ics  only  at  very  shallow  depths  (less  than  122  m  for  the  Pahute  Mesa  calculations). 
Given  that  the  tuff  model  was  chosen  to  be  a  simple  "ensemble"  average  of  structures 
from  a  number  of  shots  in  the  Pahute  Mesa  area,  this  may  reflect  the  possibility  that 
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Figure  4.17a.  Displacement  teleseismic  P  waveforms  (TOA=20*)  for  the 
three  2D  calculations  and  the  ID  RVP  convolved  with 
the  Green’s  function.  Note  the  three  2D  calculations 
predict  very  similar  P-waves  while  the  ID  synthetic 
predicts  a  large  pP  arrival.  Takeoff  angle  20* 
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Figure  4.17b.  P  waveforms  (TOA  =  20*)  from  4.17a  convolved  with 
KS36000  and  attenuation  operator  (t*  =  0.8  sec).  2D 
synthetics  all  predict  very  similar  P-waves,  while  the  ID 
synthetic  predicts  a  much  larger  “c”  phase. 
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Figure  4.18a.  Far-field  SV  waves  (TOA  =  10*)  from  the  three  2D 
calculations  compared  with  the  ID  elastic  point  source 
explosion  model.  Takeoff  angle  10*.  Note  the  direct  SV 
wave  (at  0.4  sec)  in  the  2D  calculations  that  is  absent  in 
the  ID  synthetic.  2D  synthetics  are  largely  alike  except 
for  the  amplitude  of  the  direct  SV  wave. 
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Figure  4.18b.  Far-field  SV  waves  (TOA  =  20*)  from  the  three  2D 
calculations  compared  with  the  ID  elastic  point  source 
explosion  model.  Takeoff  angle  20*.  Note  the  direct  SV 
wave  (at  0.4  sec)  in  the  2D  calculations  that  is  absent  in 
the  ID  synthetic.  2D  synthetics  are  largely  alike  except 
for  the  amplitude  of  the  direct  SV  wave. 
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Figure  4.18c.  Far-field  SV  waves  {TOA  =  30*)  from  the  three  2D 
calculations  compared  with  the  ID  elastic  point  source 
explosion  model.  Takeoff  angle  30*.  Note  the  direct  SV 
wave  (at  0.4  sec)  in  the  2D  calculations  that  is  absent  in 
the  ID  synthetic.  The  15  MPa  prestress  case  exhibits 
the  smallest  SV  amplitude,  the  no  prestress  has  the 
greatest  largest  SV  amplitude,  and  the  7.5  MPa  prestress 
and  ID  elastic  synthetics  are  intermediate. 
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Figure  4.19a.  Far-field  P  and  SV  waveforms  are  shown  to  the  same 
scale  for  comparison.  Takeoff  angle  30*.  Top  to  bottom: 
2D  no  prestress,  150  bar  prestress,  75  bar  prestress,  ID. 
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Figure  4.19b.  Linear-linear  spectra  of  waveforms  shown  in  4.19a.  Note 
that  the  P  and  SV  spectra  are  different  shapes  for  the 
2D  calculations,  unlike  the  ID  elastic  point  source  model. 
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there  is  no  "average"  structure  when  nonlinear  calculations  are  concerned.  We  suggest 
four  possible  explanations.  1)  The  material  between  122  meters  and  the  water  table  is 
weaker  than  that  used  in  the  model  such  that  spall  would  extend  to  greater  depth.  2) 
The  presence  of  weak  thin  beds  of  ash  fall  tuff,  or  other  thin  layers  of  weak  material 
may  act  as  regions  for  localized  tensile  failure.  3)  Lateral  heterogeneities  in  the 
material  above  the  shot  may  initiate  instabilities  that  localize  spall  and  produce  non- 
uniform  spalling.  4)  The  nonlinear  properties  of  the  materials  above  the  shot 
attenuated  too  much  energy  from  the  free-surface  reflection.  In  each  case,  detailed 
ground  motion  observations  and  numerical  experiments  would  help  to  test  the 
hypothesis.  In  particular,  the  detailed  modeling  of  the  inter-bedded  nature  of  the 
Pahute  Mesa  geologic  environment  would  test  whether  the  weak  beds  can  localize  ten¬ 
sile  failure  along  specific  horizons  and  extend  spall  to  greater  depth.  There  are  a  few 
limited  data  sets  of  surface  and  buried  accelerometer  records  available  to  help  control 
such  simulations. 

The  regional  Pn  observations  at  KNB  and  MNV  show  that  the  pPn  is  largely 
absent  for  Pahute  Mesa  tuff.  This  is  in  contrast  to  the  prediction  from  elastic  theory 
but  in  agreement  with  the  nonlinear  simulations.  This  indicates  that  the  simulations 
are  properly  modeling  the  nonlinear  effects  at  least  for  the  regional  Pn  waveform. 

In  the  spectral  domain,  the  synthetics  demonstrate  that  the  overburied  explosion  is 
enriched  in  high  frequencies,  and  depleted  in  low-frequencies  with  respect  to  the  nor¬ 
mal  containment  depth.  The  observations  for  the  far-field  direct  P  wave  are  consistent 
with  cube-root  depth  scaling  of  the  comer  frequency.  The  nonlinear  far-field  P-wave 
spectra  demonstrate  that  the  pP  is  not  a  spectral  replica  of  the  direct  P,  and  that  the 
nulls  predicted  by  linear  propagation  from  a  point  source  are  not  at  the  same  locations 
and  are  much  less  prominent. 

For  the  2D  axisyirunetric  simulations,  the  bulk  (but  not  all)  of  the  SV  energy 
appears  to  originate  with  the  free-surface  interact’ on  as  with  the  elastic  ID  explosion 
model.  However,  the  two-dimensional  far-field  SV  radiation  is  enhanced  in  low  fre¬ 
quencies  relative  to  the  ID  explosion  model  and  deficient  in  high  frequencies.  The 
same  high/low  frequency  effect  is  true  of  the  P-wave  radiation  but  much  less  pro¬ 
nounced.  Therefore,  we  can  not  explain  the  high/low  frequency  change  in  the  SV 
radiation  by  a  simple  comer  frequency  shift  in  an  explosive  point  source  radiator.  A 
more  complicated  model  of  the  explosion  source  is  required  to  explain  the  frequency 
content  of  the  far-field  SV. 
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V.  TWO-DIMENSIONAL  SOURCE  REGIONAL  SYNTHETICS 
5.1  Introduction 

Regional  wave-guide  phases  such  as  Pg  and  Lg  may  be  interpreted  in  three  basic 
ways;  integration  over  modes,  integration  over  generalized  rays,  or  integration  in  the 
frequency  wavenumber  domain.  In  the  modal  domain,  the  phases  are  viewed  as  a 
superposition  of  modes  with  similar  group  velocity  and  cover  a  band  of  phase  veloci¬ 
ties.  Because  the  excitation  of  each  mode  has  a  complicated  depth  dependence  and 
because  the  number  of  modes  grows  with  increasing  frequency,  we  can  not  simply 
synthesize  the  high  frequency  regional  seismogram  with  a  few  well  chosen  modes 
excited  by  a  point  source  at  a  chosen  depth.  Fimhermore,  as  is  demonstrated  by  Day, 
et  al.  (in  preparation)  the  locked  mode  approximation  to  modal  summation  synthesis 
is  shown  to  fail  for  regional  seismograms  with  attenuation.  Asymptotically,  we  can 
view  the  Pg  and  Lg  arrivals  as  bundles  of  critically  reflected  rays  making  multiple 
reflections  (rays)  in  the  crustal  waveguide.  Haskell  (1966)  showed  that  the  Pg  modes 
are  asymptotic  to  the  phase  and  group  velocities  predicted  from  such  a  ray-theoretic 
approach.  Furthermore,  Haskell  argued  that  because  of  conversion  at  interfaces  Pg 
constitutes  a  set  of  leaky  modes  that  suffers  from  loss  of  SV  energy  into  the  mantle. 
Finally,  in  the  frequency-wavenumber  domain  we  can  view  Pg  and  Lg  as  slowness 
bands  of  the  layer  response  that  constructively  interfere  to  produce  arrivals  with  a 
nearly  constant  group  velocity  as  a  function  of  distance.  Consequently,  we  can  not  syn¬ 
thesize  the  Pg  and  Lg  waveforms  with  a  single  slowness  or  take  off  angle  as  was  done 
to  produce  the  far-field  P  and  S  waveforms. 

Although  the  modal  and  ray-theoretic  methods  add  a  sense  of  intuition  to  the  pro¬ 
cess  of  generating  synthetics,  the  wavenumber  integration  synthesis  is  the  most  exact 
method  of  analysis  under  the  approximation  that  the  Earth  consists  of  laterally  homo¬ 
geneous  layers.  Because,  intuition  can  often  be  fooled  by  complex  wave  phenomena, 
we  have  used  the  more  exact  method  of  wave  synthesis  to  test  the  hypothesis  that  non¬ 
linear  2D  explosion  effects  will  produce  a  signature  in  regional  phases  and  their  spec¬ 
tra. 

The  method  used  in  the  following  section  to  generate  regional  synthetics  is  based 
upon  an  application  of  the  elastodynamic  representation  theorem  (see  Aki  and 
Richards,  1980)  and  is  exact  under  the  layered  structure  approximation.  The 
mathematical  details  are  contained  in  an  appendix.  In  essence,  the  displacements  and 
tractions  are  propagated  from  the  finite  difference  grid  to  regional  distances  using  a 
surface  boundary  integral  containing  the  wavenumber  integration  Green’s  functions. 
For  the  Pahute  Mesa  synthetics,  the  regional  velocity  model  BR2  is  used  with  the 
upper  most  layers  altered  to  exactly  replicate  the  model  used  in  the  simulations  (Table 
4.4).  In  the  case  of  the  Shagan  River  simulations,  a  structure  for  Eastern  Kazakh  from 
Stevens  (1986)  was  used  to  compute  Green’s  functions.  The  calculations  were  per¬ 
formed  on  a  CRAY-2  supercomputer.  A  typical  synthetic  at  100  km  distance,  with  a 
bandwidth  of  0  to  2  Hz,  required  approximately  90  minutes  of  CPU  time.  This  is 
largely  due  to  the  large  number  of  Green’s  functions  that  must  be  computed  at  the 
depths  and  distances  required  for  the  boundary  integral  method  of  seismogram 
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synthesis.  The  code  makes  partial  use  of  the  vector  processing  capability  of  the  super¬ 
computer.  Calculations  for  broadband  seismograms  scale  quadratically  with  distance 
and  linearly  with  bandwidth.  For  example,  a  doubling  of  range  increases  computation 
time  by  a  factor  of  four  while  doubling  the  bandwidth  simply  doubles  computation 
time.  The  nonlinear  finite  difference  calculations  are  considered  valid  up  to  5  Hz. 

5.2  Pahute  Mesa  Synthetics 

The  simulations  described  in  Section  4.3  were  propagated  as  broadband  seismo¬ 
grams  to  a  distance  of  100  km  for  frequencies  up  to  2  Hz.  Figures  5.1a  and  b  show 
regional  synthetics  for  2D  axisymmetric  calculations  compared  to  those  for  point 
source  explosion  responses  convolved  with  an  RVP.  We  see  that  the  ID  point  source 
peak  amplitude  is  larger  than  the  2D  amplitude  in  all  three  depths  of  burial.  However, 
the  general  character  of  the  seismograms  is  similar  for  the  ID  point  sources  and  the 
2D  distributed  sources.  The  2D  cratering  synthetic  has  a  more  prominent  fundamental 
Rayleigh  wave  arrival  at  55-65  seconds  than  the  ID  synthetic.  Somewhat  surprisingly, 
the  680  and  980  meter  2D  synthetics  show  more  impulsive  initial  P  waves  that  the  ID 
point  sources. 

The  surface  low  velocity  layers  in  the  model  are  responsible  for  two  very  slow 
Airy  phases  at  80-90  seconds  and  100-110  seconds.  These  slow  arrivals  show  the 
largest  differences  between  the  2D  and  ID  synthetics  and  illustrate  that  the  2D  sources 
have  nonlinear  phenomena  extending  up  to  the  near  surface  to  excite  waves  in  these 
layers.  The  low  velocity  surface  layers  responsible  for  these  Airy  phases  are  rarely 
continuous  in  the  Basin  and  Range  for  more  than  50  km  so  we  would  not  expect  these 
parts  of  the  seismogram  to  be  robust  features  of  the  propagation.  However,  the  energy 
of  these  modal  waves  must  be  scattered  into  other  portions  of  the  seismic  field  and 
contribute  energy  to  other  portions  of  the  seismogram  such  as  Lg  or  teleseismic  P 
coda. 

The  whole  record  spectra  for  the  2D  synthetics  are  shown  in  Figures  5.2a  and  b. 
The  680  and  980  spectra  are  virtually  indistinguishable  between  0.6  and  2  Hz.  Below 
0.6  Hz,  the  680  simulation  is  enhanced  in  low  firequencies  with  respect  to  the  over¬ 
buried  simulation.  The  spectra  for  the  cratering  synthetics  are  displaced  from  the  other 
spectra  for  clarity  and  are  clearly  different  from  the  contained  synthetics  above  0.5  Hz. 

Figures  5.3a  and  5.3b  show  multiple  bandpass  filtered  seismograms  of  the  2D 
vertical  simulated  seismograms  from  Figure  5.1a.  The  1  octave  bandpasses  demon¬ 
strate  that  the  680  meter  deep  source’s  surface  waves  between  0.03  and  0.5  Hz  are  70 
to  100%  larger  than  the  980  meter  synthetics.  Murphy’s  (1977)  cube-root  depth  scal¬ 
ing  would  predict  12%.  The  Pg  and  Lg  amplitudes  in  the  0.5  to  2.0  Hz  bandwidth  are 
much  the  same  amplitudes  between  the  two  depths  of  burial.  The  multiple  bandpass 
filters  therefore  show  that  the  low  frequency  enhancement  of  the  normal  buried  explo¬ 
sion  is  primarily  in  the  slow  surface  waves.  Simple  point  source  explosion  scaling  can 
not  account  for  these  excitation  differences. 

Spectra  of  the  synthetic  Pg  show  no  systematic  differences  between  the  two 
depths  of  burial,  while  a  spectrum  of  the  Lg  window  shown  in  Figure  5.4  indicates 
that  the  normal  depth  of  burial  synthetic  is  enriched  in  low  frequencies  (below  1  Hz) 
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Figure  5.1a.  Comparison  of  broadband  velocity  synthetics  for  the 
Pahute  Mesa  tuffs  calculations  at  100  km.  Vertical 
component,  0  to  2  Hz  bandwidth. 
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Figure  5.1b.  Comparison  of  broadband  velocity  synthetics  for  the 

Pahute  Mesa  tuffs  calculations  at  100  km.  Radial 
component,  0  to  2  Hz  bandwidth. 
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Figure  5.2a.  Whole  record  spectra  for  the  2D  Pahute  Mesa  tuff  synthetics  in  5.1a.  Note  the 

low-frequency  difference  between  the  680  and  980  meter  depth  synthetics. 
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Figure  5.2b.  Whole  record  spectra  for  the  2D  Pahute  Mesa  tuff  synthetics  in  5.1b.  Note  the 

low-frequency  difference  between  the  680  and  980  meter  depth  synthetics. 
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Figure  5.3a,  Bandpass  filtered  seismograms  for  the  2D  680  meter  depth  of  burial  vertical 

component  synthetic.  Traces  are  from  top  to  bottom,  broadband  velocity,  0.03-0.06 
H/,  0.06-0.12  Hz,  0.12-0.25,  and  0.25-0.5  Hz. 
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Figure  5.3b.  Bandpass  filtered  seismograms  for  the  2D  980  meter  depth  of  burial  vertical 

component  synthetic.  Traces  are  from  top  to  bottom,  broadband  velocity,  0.03-0.06 
Hz,  0.06-0.12  Hz,  0.12-0.25,  and  0.25-0.5  Hz. 
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Figure  5.4.  Whole  record  spectral  ratio  of  the  680  and  980  meter 

depth  synthetics.  The  680  meter  depth  source  is 
enhanced  in  frequencies  below  1  Hz. 
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and  depleted  in  high  frequencies  (above  1  Hz)  with  respect  to  the  overburied  synthetic. 

In  order  to  examine  this  frequency  dependence  above  2  Hz  and  to  test  the  fre¬ 
quency  dependence  upon  distance,  three  additional  narrow  bandwidths  were  examined, 
1.25-1.56,  2.50-2.81,  and  4.69-5.0  Hz  at  a  distance  of  320  km.  A  comparison  of  the 
whole  record  spectra  are  shown  in  Figure  5.5.  The  average  log-amplitudes  for  the 
three  bandwidths  show  that  for  the  two  0.31  Hz  bands  centered  on  1.41  and  2.66  Hz, 
the  two  depths  of  burial  are  indistinguishable.  However,  at  the  highest  band  centered 
at  4.85  Hz,  the  680  meter  depth  of  burial  synthetic  is  depleted  by  a  factor  of  2  with 
respect  to  the  980  meter  depth  synthetic.  The  two  spectral  averages  are  distinct  at  the 
2  a  confidence  level.  The  time  domain  seismogram  limited  to  the  4.69  to  5  Hz 
bandwidth  indicates  that  the  energy  is  predominately  P  with  high  group  velocities. 
This  result  is  qualitatively  consistent  with  the  far-field  P  spectra  of  Section  4.3.2  that 
demonstrated  that  the  far-field  P  wave  from  the  over  buried  simulation  was  enriched  in 
high  frequencies  in  accordance  with  cube-root  depth  scaling. 

In  order  to  examine  the  excitation  differences  below  2  Hz  that  can  not  be 
explained  by  cube-root  depth  scaling  of  point  explosion  sources,  we  computed  point 
explosion  (jieen’s  functions  for  a  range  of  depths.  From  symmetry  arguments,  we  can 
argue  that  the  2D  axisymmetric  source  may  to  lowest  order  be  represented  by  a  verti¬ 
cal  force,  Fz,  and  a  moment  tensor  with  three  non-zero  elements,  ,  Myy ,  and  • 
The  axisymmetric  source  has  M„=Myy.  The  explosion  part  is  given  by  M/  =  (M^  + 
Myy  +  M„  )/3  =  2/3  Mxx  +  1/3  M^^ .  The  non-isotropic  part  is  given  by  Mciyo  -  ^xx 
-  ,  and  is  termed  the  compensated  linear  vector  dipole  (CLVD).  It  is  equivalent  to 

vertical  tension  compensated  by  horizontal  compression.  Alternatively,  the  CLVD 
source  is  equivalent  to  the  superposition  of  two  45°  thrust  faults  with  fault  strikes 
oriented  at  90°  to  each  other.  It  should  be  emphasized  that  these  three  "canonical" 
sources,  Fz,  explosion,  and  CLVD,  are  only  the  first  three  terms  in  an  infinite  expan¬ 
sion  that  is  guaranteed  to  converge  only  at  low  frequencies  (wavelengths  longer  than 
the  source  dimensions). 

Conservation  of  linear  momentum  guarantees  that  Fz  has  no  net  value  (low  fre¬ 
quency).  Day,  et  al.  (1983)  have  argued  that  in  the  long  wavelength  limit,  spall  can 
be  described  by  a  point  force  model  with  with  negligible  low-frequency  amplitude. 
However,  the  vertical  point  force  and  (XVD  source  from  spall  may  radiate  short 
period  waves.  The  vertical  point  force  arises  in  the  simulations  from  the  gravitational 
attraction  that  keeps  the  grid  from  flying  off  into  space.  This  force  can  be  computed 
directly  from  the  monitored  tractions  and  displacements  on  the  cylindrical  surface  and 
is  found  to  be  small  compared  to  the  net  first  order  moments.  A  detailed  analysis  in 
terms  of  these  moments  will  be  examined  in  a  future  report,  so  we  do  not  elaborate 
further  on  the  details  of  the  expansion  except  to  show  the  Green’s  functions  for  the 
(XVD  source  as  a  function  of  depth  compared  to  the  complete  2D  axisymmetric  syn¬ 
thetics. 

Figures  5.6a,  b,  and  c  show  the  complete  2D  synthetics  (680  meter)  compared  to 
the  explosion  and  (XVD  Green’s  functions  in  three  bandwidths  of  interest  In  the  0.5 
to  1.25  Hz  bandwidth,  we  see  that  the  2D  synthetic  (at  bottom)  is  similar  to  the  explo¬ 
sion  Green’s  function  convolved  with  an  RDP  except  for  the  slow  arriving  surface 
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Figure  5.5.  Average  spectral  levels  of  synthetics  for  680  and  980 

meter  depth  sources  in  three  frequency  bands.  The  two 
depths  of  burial  are  indistinguishable  at  1.4  or  2.5  Hz. 
The  980  meter  synthetic  is  enriched  for  frequencies  above 
4.5  Hz  with  respect  to  the  680  meter  synthetic. 
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Figure  5.6a.  Bandpass  filtered  synthetics  and  Green’s  functions  in  the 

0.5  to  1.25  Hz  bandwidth.  From  bottom  to  top;  2D  680 
meter  depth  synthetic,  the  ID  synthetic  at  680  meters, 
the  explosion  Green’s  function  for  a  source  at  680  meters, 
the  CLVD  Green’s  function  for  a  source  at  680  meters, 
and  a  CLVD  Green’s  function  for  a  source  at  200  meters 
depth.  Note,  that  the  late  arriving  surface  waves  are 
best  excited  by  a  CLVD  source  and  are  not  adequately 
excited  by  an  explosive  source  at  the  working  point. 
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Figure  5.6c.  Bandpass  filtered  synthetics  and  Green’s  functions  in  the 

1.25  to  2.0  Hz  bandwidth.  From  bottom  to  top;  2D  680 
meter  depth  synthetic,  the  ID  synthetic  at  680  meters, 
the  explosion  Green’s  function  for  a  source  at  680  meters, 
the  CLVD  Green’s  function  for  a  source  at  680  meters, 
and  a  CLVD  Green’s  function  for  a  source  at  200  meters 
depth.  Note  that  neither  the  point  explosion  nor  the 
CLVD’s  excite  the  Pg  wavetrain  observed  in  the  2D 
synthetic.  The  Lg  wavetrain  is  also  longer  than  either 
the  point  explosion  or  the  CLVD  Green’s  function. 
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waves  and  some  energy  at  40  seconds.  The  CLVD  Green’s  functions  for  depths  of 
680  and  200  meters  are  rich  in  energy  at  these  times  and  weak  at  other  arrival  times. 
Clearly  the  CLVD  source  is  an  efficient  low  phase  velocity  SV  radiator  without  pro¬ 
ducing  substantial  Pg  waves  or  the  surface  waves  that  arrive  at  45  to  SO  seconds.  In 
the  0.75  to  1.5  Hz  and  1.25  to  2.0  Hz  bandwidths  we  see  that  the  explosion  Green’s 
function  can  account  for  the  majority  of  the  2D  synthetic  except  for  the  surface  wave 
that  arrives  at  80  seconds  and  some  of  the  high  frequency  impulsive  arrivals  in  the  Pg 
coda  (20  to  40  seconds).  The  CLVD  source  is  an  efficient  radiator  for  the  surface 
wave  arriving  at  80  seconds.  A  shallow  CLVD  source  (2(X)  meters)  is  an  even  better 
radiator  than  a  deeper  source  (680  meters). 

Consequently,  we  can  see  that  the  2D  synthetics  differ  from  a  simple  explosion 
point  source.  It  is  not  possible  to  model  the  2D  explosion  source  with  a  single  explo¬ 
sion  Green’s  function  convolved  with  an  RDP.  The  relative  excitation  of  Pg,  Lg,  and 
surface  modal  waves  as  a  function  of  frequency  demonstrate  that  a  higher  order 
description  of  the  source  beyond  a  point  explosion  is  necessary.  The  normal  buried 
source  is  enriched  in  shallow  modal  excitation  at  low  frequencies  (below  12  Hz)  rela¬ 
tive  to  either  a  point  source  or  the  overburied  source.  The  deeper  source  is  enriched  in 
high  frequencies  in  the  Lg  (above  1  Hz)  and  Pg  (above  4  Hz).  The  enrichment  of  high 
frequency  P  waves  for  overburial  is  consistent  with  cube-root  depth  scaling. 

The  very  simple  examination  of  the  CLVD  Green’s  functions  compared  to  the  2D 
synthetics  indicates  that  a  superposition  of  CLVD  and  explosion  sources  may  model 
much  of  the  observed  excitation  behavior  of  the  2D  synthetics.  The  CLVD  source 
would  have  a  different  spectral  content  (and  therefore  time  function)  that  provides  the 
enhancement  of  the  shallow  source  with  respect  to  the  deeper  source  in  the  0.5  to  1.5 
Hz  bandwidth.  An  obvious  physical  mechanism  for  this  "signal"  is  the  spall  behavior 
of  the  sources  as  a  function  of  scaled  depth.  The  deeper  source  has  less  spall;  the 
spall  signal  is  weaker  and  shorter  in  duration.  Asymmetrical  expansion  of  the  source 
such  as  less  (or  more)  expansion  in  the  horizontal  direction  than  the  vertical  direction 
could  also  result  in  a  net  CLVD  component  to  the  source. 

As  was  pointed  out  earlier,  the  Pahute  Mesa  simulations  exhibited  spall,  but  it 
was  confined  to  a  weak  shallow  surface  layer.  Observations  indicate  that  spall  typi¬ 
cally  occurs  at  greater  depths  than  in  the  simulations,  therefore,  we  would  expect  that 
the  spall  signal  would  be  larger  in  reality  than  was  predicted  by  the  2D  calculations. 
The  CLVD  source  is  an  efficient  radiator  of  P  waves  at  steep  angles  of  incidence  and  a 
weak  P-wave  radiator  at  shallow  angles  of  incidence.  Therefore,  a  larger  CLVD  may 
be  the  source  component  needed  to  improve  the  teleseismic  P  waveform  simulations 
and  provide  a  mechanism  that  can  explain  the  high/low  frequency  observations  of 
Murphy  and  Bennett  (1982)  and  Taylor,  et  al.  (1988). 

5.3  Shagan  River  Synthetics 

The  model  used  to  compute  regional  synthetics  for  the  Shagan  River  calculations 
is  listed  in  Table  5.3,  and  is  taken  from  Stevens  (1986)  with  the  minor  addition  of  a  2 
km  thick  layer  identical  to  that  used  in  the  finite  difference  calculations.  For  com¬ 
parison,  Figure  5.7  compares  the  Eastern  Kazakhstan  P-velocity  models  of  Leith 
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(1987),  and  Priestley,  et  al.  (1988)  with  the  Stevens  (1986)  model.  These  models  are 
characterized  by  a  thick  high  velocity  crust  and  a  high  velocity  sub-Moho  velocity. 
The  Priestley,  et  al.  crustal  models,  based  on  vertical-to-radial  teleseismic  P-wave 
transfer  functions,  are  characterized  by  nearly  continuous  crustal  velocity  gradients. 
Leith’s  model,  based  on  Soviet  published  regional  refraction  profiles,  is  similar  to  the 
Priestley,  et  al.  models  with  no  prominent  "Conrad"  discontinuity.  Stevens’  model 
based  on  inversion  of  explosion  generated  waves  recorded  at  MAIO,  and  other  pub¬ 
lished  Soviet  profiles  places  most  of  the  crustal  velocity  increase  in  a  mid-crustal  tran¬ 
sition  zone. 


TABLE  5.3  Eastern  Kazakhstan  Structure 


Modified  from  Stevens  (1986) 


h(km) 

a(km/s) 

3(km/s) 

2. 

5.02 

2.79 

1. 

5.4 

3.0 

2.488 

5.9 

3.3 

10.976 

6.1 

3.4 

5.488 

6.308 

3.541 

5.488 

6.597 

3.703 

5.488 

6.736 

3.781 

5.564 

6.782 

3.807 

6.504 

6.795 

3.814 

8.006 

8.147 

4.573 

9.359 

8.138 

4.568 

10.94 

8.106 

4.55 

12.78 

8.065 

4.527 

14.95 

8.047 

4.517 

17.47 

8.07 

4.53 

20.42 

8.117 

4.556 

23.88 

8.154 

4.577 

27.91 

8.161 

4.581 

32.63 

8.145 

4.572 

38.14 

8.12 

4.558 

oo 

8.101 

4.547 

p(gm/cc) 

(2u 

depth 

2.7 

100. 

0 

2.7 

150. 

2 

2.7 

200. 

3 

2.7 

600. 

5.488 

2.702 

525.9 

16.464 

2.807 

500. 

21.952 

2.858 

450. 

27.44 

2.875 

400. 

32.928 

2.879 

350. 

38.492 

3.372 

179.5 

44.996 

3.369 

167.5 

53.002 

3.358 

159.4 

62.361 

3.343 

153.8 

73.301 

3.336 

150.3 

86.081 

3.345 

148.7 

101.031 

3.361 

148.1 

118.501 

3.375 

147.9 

138.921 

3.378 

147.4 

162.801 

3.372 

146.8 

190.711 

3.363 

146.2 

223.341 

3.356 

145.8 

261.481 

The  synthetics  for  an  epicentral  distance  of  100  km  are  shown  in  Figures  5.8a  and 
b.  As  with  the  Puhute  Mesa  synthetics  of  the  previous  section,  each  2D  calculation  is 
paired  with  a  ID  calculation.  The  largest  amplitude  phase  on  the  broadband  as  well  as 
bandpass  filtered  seismograms  (Figures  5.9a,  b,  and  c)  up  to  about  1.2- 1.4  Hz,  is  a  fun- 
daniental  mode  (LR  long-period  or  Rg  short-period)  Airy  phases  with  group  velocities 
between  2.7  and  3.0  km/s.  The  impulsive  direct  crustal  P  (Pg)  wave  is  followed  by 
SV  arrivals  with  "group"  velocities  between  3.5  and  2.7  km/s.  Maximum  amplitudes 
from  the  bandpass  filtered  seismograms  of  Figures  5.9  and  Fourier  spectra  are  com¬ 
pared  in  Figures  5.10a,  b,  and  c. 
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Figure  5.8a.  Broadband  velocity  synthetics  for  the  Shagan  River 

calculations  at  200,  680,  and  980  meter  depths.  2D  and 
ID  synthetics  are  displayed  together  for  comparison. 
Vertical  component,  100  km,  0  to  2  Hz  bandwidth. 
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Figure  5.8b.  Broadband  velocity  synthetics  for  the  Shagan  River 
calculations  at  200,  680,  and  980  meter  depths.  2D  and 
ID  synthetics  are  displayed  together  for  comparison. 
Radial  component,  100  km,  0  to  2  Hz  bandwidth. 
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Figure  5.9a.  Bandpass  filtered  Shagan  River  2D  synthetics  (left)  and 

ID  synthetics  (rightj  for  depth  of  200  meters. 
Bandwidths  are  indicated  to  the  left. 


Figure  5.9b.  Bandpass  filtered  Shagan  River  2D  synthetics  (left)  and 
ID  synthetics  (right)  for  depth  of  680  meters. 
Bandwidths  are  indicated  to  the  left. 


Figure  5.9c.  Bandpass  filtered  Shagan  River  2D  synthetics  (left)  and 
ID  synthetics  (right)  for  depth  of  980  meters. 
Bandwidths  are  indicated  to  the  left. 
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The  whole  record  spectra  largely  reflect  the  time-domain  peak  amplitudes  and  are 
dominated  by  the  fundamental  mode  Rayleigh  up  to  1.0  Hz  for  the  2D  or  1.2  Hz  for 
the  ID  synthetics  at  depths  of  680  and  980  meters.  The  200  meter  depth  synthetics 
are  dominated  by  the  fundamental  up  to  1.5  Hz  for  the  2D  synthetic  and  up  to  the  2 
Hz  Nyquist  frequency  for  the  ID  synthetic.  The  large  amplitude  of  the  ID  200  meter 
depth  synthetic  is  a  product  of  a  large  ID  moment  (a  factor  of  2  over  the  980  meter 
synthetic)  and  strong  depth  dependence  for  the  excitation  function  of  the  fundamental 
mode  elastic  Green’s  function.  Clearly,  the  ID  synthetic  demonstrates  the  short  com¬ 
ings  of  the  point  explosion  source  model  for  the  cratering  explosion  in  this  structure. 
The  differences  between  the  2D  and  ID  contained  synthetics  (680  and  980  meters)  are 
more  subtle.  The  prominent  "notch"  in  the  whole  record  spectra  about  1.0  Hz  in  the 
2D  synthetics  is  at  1.4  Hz  in  the  ID  synthetics.  This  notch  occurs  where  the  funda¬ 
mental  mode  amplitude  decreases  as  a  function  of  increasing  frequency  and  declines 
below  the  level  of  the  higher  modes  and  SV  body  waves  which  are  increasing  with 
frequency.  Obviously  this  "notch"  will  be  a  complicated  function  of  burial  depth, 
event  size  (hence  comer  frequency),  and  distance  (hence  differential  attenuation  of 
modes  and  body  waves).  The  peak  amplitudes  taken  from  the  whole  record  do  not 
show  this  spectral  "notch"  nearly  as  well  as  the  power  spectrum  because  the  signal 
power  spectra  is  dispersed  in  the  time  domain  as  several  arrivals. 

The  P  and  SV  window  spectra  and  amplitudes  are  shown  in  Figures  5.10b  and  c. 
The  2D  P-wave  spectra  and  time  domain  peak  amplitudes  illustrate  that  the  680  meter 
depth  synthetics  have  an  enhancement  (about  a  factor  of  2)  in  the  0.4  to  0.8  Hz 
bandwidth  with  respect  to  the  980  meter  depth  synthetics.  This  is  not  the  case  for  the 
ID  elastic  synthetics  which  show  a  nearly  constant  amplitude  offset  between  the  two 
source  depths  of  about  40%  up  to  1.6  Hz.  The  ID  spectra  simply  reflect  the  cube-root 
depth  scaling  of  the  ID  RVP’s.  This  is  clearer  when  we  plot  the  2D  spectra  refer¬ 
enced  to  the  ID  spectra  in  Figure  5.11.  Only  the  two  depths,  680  and  980  meters,  are 
shown.  First,  note  that  the  P-wave  low  frequencies  for  the  680  meter  depth  event  are 
enhanced  with  respect  to  the  980  meter  depth  event.  Both  contained  P-waves  show  a 
tendency  for  there  to  be  a  loss  of  high  frequencies  above  1.6  Hz.  Second,  note  that 
the  680  meter  depth  event  whole  record  spectra  show  an  enhancement  in  the  0.2  to  0.8 
Hz  frequency  bandwidth.  There  is  little  systematic  difference  between  the  two  whole 
record  spectra  above  1.0  Hz.  Third,  the  SV  wavetrain  (2.7-3.5  km/s)  spectra  show  lit¬ 
tle  difference  between  the  two  depths  of  burial.  The  2D  synthetics  are  enhanced  in  the 
1.2  to  1.6  Hz  frequency  bandwidth  with  respect  to  the  ID  explosion  source  model. 

As  with  the  Pahute  Mesa  synthetics,  we  present  in  Figures  5.12a,  b,  and  c,  a  com¬ 
parison  at  three  bandwidths  of  the  2D  synthetic  with  the  explosion  and  CLVD  Green’s 
functions.  In  the  0.5  to  1.25  Hz  bandwidth,  there  is  an  apparent  delay  in  the  funda¬ 
mental  Rayleigh  wave  of  the  2D  source  with  respect  to  the  explosion  or  CLVD 
Green’s  functions.  Actually  this  represents  the  lower  frequency  nature  of  the  2D  fun¬ 
damental  relative  to  the  ID  point  source  model.  There  can  not  be  much  CLVD  in  the 
2D  synihetic  in  this  frequency  band.  In  the  0.75  to  1.5  Hz  bandwidth,  the  2D  syn¬ 
thetic  has  arriving  energy  in  the  30  to  40  second  time  window  that  is  absent  in  the 
explosion  Green’s  function  and  easily  accounted  for  by  a  CLVD  source.  In  the  1.25  to 
2.0  Hz  bandwidth,  the  arrivals  in  the  30  to  40  second  window  can  easily  be  explained 
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SHAGAN  RIVER  REGIONAL  SYNTHETICS  100km  0.5-1.25Hz 
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Figure  5.12a.  Shagan  River  680  meter  depth  of  burial  2D  synthetic  and 
several  Green’s  functions  in  the  0.5-1.25  bandwidth. 
From  bottom  to  top  are:  the  680  meter  depth  2D 
synthetic,  the  680  meter  depth  explosion  Green’s  function 
convolved  with  the  ID  RDP,  the  680  meter  depth 
explosion  Green’s  function,  the  CLVD  Green’s  function  at 
680  meters  depth,  and  the  CLVD  Green’s  function  at  200 
meters. 
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SHAGAN  RIVER  REGIONAL  SYNTHETICS  100km  0.75-1.5Hz 


Figure  5.12b.  Shagan  River  680  meter  depth  of  burial  2D  synthetic  and 
several  Green’s  functions  in  the  0.75-1.5  bandwidth. 
From  bottom  to  top  are:  the  680  meter  depth  2D 
synthetic,  the  680  meter  depth  explosion  Green’s  function 
convolved  with  the  ID  RDP,  the  680  meter  depth 
explosion  Green’s  function,  the  CLVD  Green’s  function  at 
680  meters  depth,  and  the  CLVD  Green’s  function  at  200 
meters.  Note  the  SV-Lg  arrival  on  the  2D  trace  that 
resembles  the  CLVD  Green’s  function. 
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Figure  5.12c. 


Shagan  River  680  meter  depth  of  burial  2D  synthetic  and 
several  Green’s  functions  in  the  1.25-2.0  bandwidth. 
From  bottom  to  top  are;  the  680  meter  depth  2D 
synthetic,  the  680  meter  depth  explosion  Green’s  function 
convolved  with  the  ID  RDP,  the  680  meter  depth 
explosion  Green’s  function,  the  CLVD  Green’s  function  at 
680  meters  depth,  and  the  CLVD  Green’s  function  at  200 
meters.  Note  the  SV-Lg  arrival  on  the  2D  trace  that 
resembles  the  CLVD  Green’s  function. 
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with  a  Q.VD  source. 

To  satisfy  these  observations,  the  CLVD  source  must  be  peaked  in  the  1  to  2  Hz 
bandwidth  to  explain  the  2D  Shagan  River  synthetics  rather  than  in  the  0.5  to  1  Hz 
bandwidth  as  with  the  Pahute  Mesa  synthetics.  However,  the  situation  is  not  as  simple 
as  with  the  Pahute  Mesa  synthetics.  There  is  a  need  to  either  reduce  the  0.5  to  1  Hz 
explosion  energy  relative  to  the  point  source  or  to  provide  some  sort  of  interference 
between  the  point  explosion  and  the  CLVD  (or  higher  order  moments)  in  this  fre¬ 
quency  band.  The  fundamental  energy  arriving  at  40  seconds  is  obviously  missing  in 
the  2D  synthetic  relative  to  the  point  explosion  or  the  CLVD  Green’s  functions. 

In  summation,  we  see  that  the  Shagan  River  simulations  indicate  that  even  a  sim¬ 
ple  source  structure  ma\  yield  a  2D  seismogram  that  can  not  be  properly  modeled  by  a 
point  explosion.  The  Pg,  Lg,  and  fundamental  all  show  evidence  of  a  complex  source. 
The  Lg  amplitudes  are  enhanced  relative  to  the  point  explosion  source,  and  the  Pg  and 
fundamental  waveforms  show  a  frequency  dependence  that  suggests  interference 
between  multiple  sources. 

5.4  Prestress  -  PILEDRIVER  Synthetics 

The  granite  calculations  presented  here,  with  and  without  prestress,  have  been 
previously  discussed  by  Day,  et  al.  (1986).  They  synthesized  long  period  teleseismic 
Rayleigh  waves  from  the  results  of  the  finite  difference  equations  and  demonstrated 
that  prestress  could  effectively  reduce  or  reverse  Rayleigh  waves  from  an  explosion. 
These  finite  difference  calculations  were  in  effect  simulations  of  the  PILEDRIVER 
event  (61  kt)  located  in  Climax  Stock  (depth  of  463  m)  at  NTS.  Near-field  motion 
was  used  to  constrain  the  2D  calculations  for  the  non  prestress  simulations  and  the 
results  are  a  set  of  well  constrained  calculations  that  properly  model  spall  as  well  as 
the  observed  near-field  RDP  estimates. 

In  addition  to  the  PILEDRIVER  calculation,  two  calculations  were  performed 
corresponding  to  horizontal  non  hydrostatic  prestresses  of  7.5  and  15  MPa  (75  and  150 
bars).  These  prestress  conditions  correspond  to  a  diagonal  stress  tensor  with  =  (5yy 
=  7.5  or  15  MPa,  =  0,  superimposed  on  the  hydrostatic  pressure  field.  Subsequent 
tectonic  stress  release  due  to  relaxation  in  the  nonlinear  zone  of  the  explosion  results 
in  an  "explosion"  source  plus  a  "tectonic"  source.  In  the  long  period  limit,  this  com¬ 
posite  source  is  composed  of  a  point  explosion  with  moment  tensor  proportional  to  the 
identity  tensor,  and  a  tectonic  source  with  moment  tensor  proportional  to  a  CLVD. 
The  two  sources  have  different  time  functions  and  therefore  introduce  an  interference 
pattern  that  is  manifested  as  reversed  and  phase  shifted  long  period  Rayleigh  waves. 

The  near  field  effects  of  prestress  are  readily  seen  in  RDP’s  estimated  from  the 
monitored  displacements  of  the  finite  difference  calculations.  Figure  5.13  shows  the 
estimated  RDP’s  from  the  radial  displacement  synthetic  record  at  a  range  of  1090 
meters  at  the  working  point  level.  The  zero  prestress  RDP  has  a  peak  near  50,(X)0 
and  a  static  value  near  20,000  {2*  10*^  Nt-M).  The  7.5  MPa  prestress  case  shows 

that  the  static  explosion  "estimated  RDP"  at  this  location  is  very  nearly  canceled  by 
the  tectonic  stress  release.  The  15  MPa  prestress  case  shows  a  complete  reversal  of 
the  "estimated  RDP". 
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Figure  5.13.  RDP  estimates  from  the  radial  displacements  of  the  finite 

difference  calculations  at  the  working  point  level.  No 
prestress  (bottom)  shows  a  typical  granite  RDP  with 
static  level  near  20,000  m"^.  7.5  MPa  prestress  (top) 
shows  a  static  level  near  zero,  while  the  15  MPa  prestress 
(middle)  shows  a  complete  reversal. 
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In  contrast  to  the  long  period  teleseismic  seismogram  the  short  period  regional 
seismogram  simulations  show  little  effects  of  prestress.  Synthetics  for  the  three  cases 
are  shown  in  Figures  5.14a,  b,  and  c  for  three  bandwidths,  0.5  to  1.25,  0.75  to  1.5  and 
1.25  to  2  Hz,  respectively.  Along  with  these  synthetics  are  Green’s  functions  for  the 
explosion  and  the  Q.VD  source  for  comparison.  The  structure  used  to  compute  the 
synthetics  is  the  BR2  model  with  the  Climax  Stock  structure  replacing  the  upper  3.95 
km  with  the  granite  structure  used  in  the  finite  difference  calculations.. 

The  three  2D  synthetics  are  remakably  like  the  explosion  Green’s  function  in  the 
short  period  band.  The  relative  amplitudes  of  the  Pg,  Lg,  and  Rayleigh  waves  are 
quite  similar  to  a  point  explosion  model.  Unlike  the  Pahute  Mesa  and  Shagan  River 
simulations,  the  2D  granite  calculations  are  modeled  fairly  well  by  the  point  explosion 
and  require  little  or  no  CLVD  source  to  match  the  Pg/Lg  ratio. 

In  order  to  demonstrate  that  the  prestress  does  produce  effects  at  longer  periods, 
we  show  in  Figure  5.15,  the  fundamental  mode  excited  by  the  three  sources.  Note  that 
periods  longer  than  10  seconds  are  reversed.  There  is  significant  miss-match  between 
the  three  seismograms  in  the  10  to  2  sec  bandwidth,  but  that  the  seismograms  become 
more  similar  with  increasing  frequency.  In  the  1  to  2  Hz  band  the  regional  seismo¬ 
gram  is  dominated  by  the  explosion  and  shows  little  evidence  of  the  tectoiuc  release. 

Modified  BR2  Structure 
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Figure  5.14a.  2D  PILEDRIVER  synthetics  and  Green’s  functions  in  the  0.5-1.25  Hz  band. 

Bottom  to  top:  the  explosion  Green’s  function,  no  prestress,  7.5  MPa  prestress,  15 
MPa  prestress,  and  the  CLVD  Green’s  function. 
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Figure  5.14c.  2D  PILEDRIVER  synthetics  and  Green’s  functions  in  the  1.25-2.0  Hz  band. 

Bottom  to  top:  the  explosion  Green’s  function,  no  prestress,  7.5  MPa  prestress,  15 
MPa  prestress,  and  the  CLVD  Green’s  function. 
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PILEDRIVER  SIMULATIONS,  FUNDAMENTAL  MODE  EXCITATION 
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Figure  5.15.  Fundamental  mode  excitation  for  the  three  2D 
PILEDRIVER  simulations.  Vertical  component,  100  km, 
0  to  2  Hz  bandwidth.  Note  that  although  the  10  second 
surface  waves  are  reversed  for  the  prestressed  cases,  the 
Rayleigh  waves  are  progressively  alike  with  increasing 
frequency  (increasing  group  velocity). 
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5.5  Conclusions 

We  see  from  the  three  sets  of  simulations,  Pahute  Mesa  tuff,  Shagan  River,  and 
prestressed  granite,  that  structure,  medium,  depth  of  burial,  and  prestress  all  have  dis¬ 
tinctly  different  effects  on  the  regional  seismogram.  However,  in  all  of  the  2D  syn¬ 
thetics  except  for  granite,  the  CLVD  component  of  the  source  is  a  significant  contribu¬ 
tor  to  SV  radiation.  This  secondary  source  modulates  P  radiation,  and  interferes  with 
the  fundamental  mode  in  the  0.25  to  2  Hz  band.  We  note  that  for  each  of  the  three 
sets  of  synthetics,  the  effects  occur  at  different  frequencies  and  are  apparent  on 
different  phases.  A  conclusion  common  to  all  the  synthetic  simulations,  except  gran¬ 
ite,  is  that  a  point  explosive  source  is  not  an  adequate  model  for  regional  explosion 
seismograms  in  the  0.25  to  2  Hz  bandwidth. 

We  saw  in  Section  III  that  he  ID  point  source  explosion  model  does  not  ade¬ 
quately  excite  regional  Lg  (Figures  3.3a,  b,  and  c).  Both  the  Pahute  Mesa  and  Shagan 
River  2D  synthetics  show  enhanced  Lg  excitation  over  the  ID  synthetics  It  appears 
that  to  first  order  introduction  of  a  CLVD  term  may  explain  some  of  the  2D  Lg  excita¬ 
tion.  However,  it  is  clear  that  the  amount  of  Lg  excitation  in  the  0.5  -  2  Hz  bandwidth 
even  for  the  2D  synthetics  falls  far  short  of  the  observations.  It  is  possible  that 
increased  spall  in  the  2D  simulations  could  contribute  to  an  increase  in  the  CLVD 
component  of  the  short  period  source  and  hence  increase  Lg  excitation  without  increas¬ 
ing  the  Pn  and  Pg  excitation. 

As  for  the  simulations  of  explosions  in  prestressed  granite,  we  see  that  the  "tec¬ 
tonic  release"  has  negligible  short  period  signature.  This  is  in  agreement  with  the  lack 
of  correlation  between  short  period  (Lg)  and  long  period  (LR)  observations,  where  the 
long  period  observations  clearly  indicate  that  tectonic  release  can  significantly  effect 
the  amplitude  and  phase  of  fundamental  Rayleigh  waves.  The  2D  simulations  suggest 
that  the  coherent  tectonic  release  has  a  gradual  onset  longer  than  two  seconds  such  that 
waves  with  period  shorter  than  two  seconds  are  virtually  unaffected. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 
6.1  Primary  Conclusions 

The  main  thrust  of  the  effort  has  been  to  examine  the  nature  of  explosion  depth 
dependence  on  the  excitation  of  regional  phases.  Several  observations  have  been  made 
over  the  last  few  years  beginning  with  Murphy  and  Bennett  (1982)  that  regional  phases 
from  explosions  at  NTS  are  lacking  in  high  frequencies  relative  to  earthquakes.  We 
have  examined  several  possible  causes  of  such  a  phenomenon.  First,  we  have  exam¬ 
ined  the  depth  dependence  of  explosion  and  earthquake  Green’s  functions.  This  work 
follows  in  the  wake  of  work  by  Bennett,  et  al.  (1987)  and  uses  a  more  accurate  tech¬ 
nique  to  model  regional  propagation,  namely  wavenumber  integration  as  opposed  to 
modal  summation.  We  find  that  although  there  are  depth  dependent  trends  in  the 
Green’s  functions  for  explosions  and  earthquakes,  ID  layered  propagation  with  simple 
point  explosion  sources  can  not  explain  the  observations. 

Point  source  explosions  in  simple  ID  layered  crustal  models  can  not  excite  ade¬ 
quate  amplitudes  of  Lg  to  match  observations.  The  nonlinear  2D  axisymmetric  finite 
difference  source  simulations  that  have  been  examined  begin  to  remedy  this  deficiency. 
Although  it  is  clear  from  the  far-field  SV  waveforms  that  the  2D  source  does  not 
directly  radiate  more  energy  than  the  ID  source.  Bennett,  et  al.  (1987)  suggested  that 
with  introduction  of  near- source  layering,  P-SV  conversion  could  be  increased  and 
increase  Lg  amplitudes.  However,  even  with  substantial  layering  near  the  source,  a 
point  explosion  source  does  not  generate  sufficient  SV  energy  in  a  layered  structure 
within  the  necessary  phase  velocity  range  of  3  to  4.0  km/s.  When  intrinsic  attenuation 
is  added  to  the  crust  (in  order  to  model  the  spatial  decay  of  Lg  at  1  Hz)  the  explosion 
point  source  synthetic  Lg  amplitude  is  practically  non-existent  beyond  300  km.  The 
"two-dimensional"  axisymmetric  nonlinear  finite  difference  simulations  examined  in 
this  study  also  suffer  from  this  deficiency  at  frequencies  below  1  Hz.  Above  1  Hz,  we 
found  that  2D  synthetics  exhibited  relative  Lg  excitation  that  required  a  component  of 
(XVD  to  the  source.  This  part  of  the  source  does  not  have  the  same  time  function  as 
the  explosion,  and  if  it  is  caused  by  spall  then  it  has  no  low  frequency  moment. 
Therefore,  we  should  expect  the  CLVD  component  of  the  source  to  be  peaked  with 
respect  to  the  explosion  source.  The  dominant  frequency  of  the  (XVD  component 
will  be  determined  by  the  time  history  of  the  nonlinear  free  surface  interaction  includ¬ 
ing  spall.  Although  the  2D  simulations  excite  more  Lg  than  the  ID  point  source,  they 
still  are  not  adequate  radiators  of  regional  SV  energy  to  agree  with  observations.  The 
regional  Green’s  functions  suggest  that  processes  that  enhance  the  CLVD  component 
of  the  source  will  generally  incrase  the  Lg/Pn  and  Lg/Pg  ratios.  These  physical 
processes  would  include,  but  not  be  limited  to,  spall  closure,  bed  separation  and  aniso¬ 
tropic  cavity  formation. 

The  Pahute  Mesa  simulations  show  that  the  primary  depth  of  burial  effects  on  the 
regional  seismogram  excitations  are  manifested  in  the  shallow  modal  content  of  the 
seismogram.  These  shallow  Rayleigh  waves  (Rg)  cannot  be  expected  to  propagate  to 
large  distances  in  a  region  such  as  the  Basin  and  Range,  but  must  be  scattered  into 
other  phases  (perhaps  Lg  and  P-coda).  The  2D  simulations  show  that  the  depth  of 
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burial  strongly  affects  the  excitation  of  these  shallow  modes  over  and  above  that 
predicted  by  a  point  explosive  source  at  the  working  point  depth.  Presumably,  the 
nonlinear  effects  of  the  source  extend  to  shallow  depths  and  give  the  event  the  appear¬ 
ance  of  being  shallower  than  the  working  point.  The  overburied  simulations  showed 
less  of  this  trend  and  therefore  we  expect  that  this  mechanism  would  contribute  to  the 
enhancement  of  low  frequencies  of  normal  burial  with  respect  to  over  burial. 

Clearly,  the  short  period  excitation  of  Lg  by  explosions  requires  some 
mechanism(s)  that  has(have)  not  been  properly  modeled.  Alternative  mechanisms 
include  scattering,  tectonic  release,  and  non-isotropic  nonlinear  effects.  One  promising 
linear  scattering  mechanism  is  the  conversion  of  short-period  Rayleigh  waves  into  Lg 
by  shallow  heterogeneity  and  topography  (McLaughlin  and  Jih,  1987).  Detailed  calcu¬ 
lations  to  verify  the  efficacy  of  this  model  have  not  yet  been  carried  out.  Although 
long-period  tectonic  release  has  not  been  shown  to  correlate  well  with  short-period  sig¬ 
nals,  there  are  observations  of  well  developed  short-period  SV  and  SH  waves  in  the  3 
to  10  km  range  (Stump  and  Johnson,  1984;  McLaughlin,  et  al.  1983).  These  shear 
waves  may  be  produced  by  some  incoherent  release  of  tectonic  strain  or  they  may  be 
due  to  some  nonlinear  interaction  of  the  explosion  with  near-source  heterogeneity. 

The  axisymmetric  finite  difference  simulations  nxxlel  nonlinear  effects  of  pore 
crush-up,  cracking,  and  spall  assuming  a  simple  layered  structure.  The  simulations 
predict  observable  direct  far-field  SV  from  the  explosion  as  well  as  from  conversion  at 
interfaces  and  spall  formation.  Simulations  have  in  the  past  stressed  simple  layered 
models  of  the  near-source  in  order  to  simplify  the  interpretation  of  the  physics 
involved.  Introduction  of  additional  layering  may  produce  more  spall  above  the  explo¬ 
sion.  Furthermore,  addition  of  lateral  heterogeneity  around  the  explosion  may  intro¬ 
duce  non-uniform  spall  and  generate  SV  and  SH  motion  not  predicted  by  the  laterally 
homogeneous  layering.  Such  a  mechanism  is  nonlinear  and  can  only  be  explored 
using  more  elaborate  nonlinear  models  with  presently  available  nonlinear  codes. 

The  CLVD  component  to  the  2D  synthetics  is  the  principal  difference  between 
the  2D  and  the  ID  models.  This  component  of  the  source  is  largely  responsible  for 
the  differences  in  regional  synthetics  between  the  2D  fiill-wave  synthetics  and  the  sim¬ 
ple  point  explosion  model.  Since  this  source  component  radiates  P  energy  largely  in 
the  vertical  direction  it  could  be  responsible  for  increasing  the  apparent  "pP"  in  telese- 
ismic  signals  while  contributing  little  "pP"  to  regional  Pn.  Furthermore,  since  it  radi¬ 
ates  SV  primarily  in  the  horizontal  direction,  this  source  component  can  be  responsible 
for  excitation  of  regional  Lg. 

Several  lines  of  evidence  suggest  that  the  Pahute  Mesa  tuff  simulations  underes¬ 
timated  the  amount  of  spall  observed  for  "typical"  events  of  their  size  and  scaled 
depth.  Spall  has  been  observed  to  occur  to  1/3  the  depth  of  burial  for  events  in  the  20 
to  150  kt  range  with  scaled  depths  of  burial  between  120  and  180  m*Y*(-l/3)  (Howard 
Patton,  personal  communication).  Spall  only  occurred  in  the  upper  most  layer  (122  m) 
for  tuff  simulations  examined  in  this  report.  Presumably,  the  low  tensile  yield  strength 
of  the  upper  most  weak  layer  served  to  confine  the  spall  to  such  a  shallow  layer. 
Analysis  of  the  Pahute  Mesa  tuff  simulations  indicates  that  the  explosion  moments  (or 
RDP’s)  are  somewhat  small  compared  to  observational  estimates  for  Pahute  Mesa 
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events.  If  the  yield  strength  of  the  tuffs  were  to  be  reduced  above  the  shots,  the 
explosion  moments  would  increase,  and  amplitudes  above  the  shots  would  increase  as 
well  as  encourage  spall  at  greater  depth.  Increased  spall  would  increase  the  CLVD 
component  of  the  source  and  therefore  increase  the  amount  of  SV-Lg  excitation. 
Increased  spall  may  contribute  to  the  "pP"  signal  at  steep  take  off  angles  and  better 
simulate  the  teleseismic  P  waveform.  At  the  same  time,  the  CLVD  component  to 
"pPn"  would  not  be  greatly  affected,  and  the  Pn  waveform  would  be  modeled  as 
before. 

6.2  Ancillary  Conclusions 

From  our  examination  of  synthetic  nwthods,  we  have  found  that  "full-wave" 
wavenumber  integration  techniques  are  required  to  properly  model  the  entire  seismo¬ 
gram  at  regional  distances  in  ±e  western  United  States.  Modal  summation  methods 
can  lead  to  erroneous  results  if  applied  for  high  frequency  synthetics  in  the  presence  of 
attenuation. 

A  crustal  shear  Q  that  increases  with  increasing  frequency  is  required  in  order  to 
properly  model  the  spectral  and  temporal  character  of  Lg  as  a  function  of  distance  in 
the  basin  and  range.  The  relative  amplitudes  of  Lg  for  synthetic  earthquake  and  explo¬ 
sion  sources  at  3  Hz  can  not  be  properly  simulated  with  a  frequency  independent  cru¬ 
stal  Q. 

Lg  synthetics  indicate  that  the  spatial  Lg  attenuation  corrected  for  Aiiy-phase  and 
cylindrical  geometrical  spreading  should  be  an  excellent  measure  of  the  "average" 
shear  Q  of  the  crust.  Similar  results  are  found  for  synthetic  Pg  attenuation  and 
geometrical  spreading.  For  the  model  examined,  the  leakage  of  Pg  energy  was 
minimal  and  the  synthetic  Pg  Q  closely  approximated  the  average  crustal  compres- 
sional  Q. 

The  frequency  dependence  of  Pn  attenuation  is  found  to  be  sensitive  to  the  gra¬ 
dient  in  the  upper  mantle  lid  just  below  the  Moho.  A  slight  velocity  gradient  below 
the  Moho  will  introduce  an  apparent  Pn  Q  that  increases  with  increasing  frequency 
even  though  a  constant  Q  intrinsic  attenuation  is  assumed. 

6.3  Recommendations 

The  results  of  this  work  suggest  that  nonlinear  scaled  depth  of  burial  effects 
should  be  observable  in  regional  seismograms.  However,  the  nonlinear  calculations 
analyzed  in  this  report  only  cover  a  limited  range  of  scaled  depths  and  yields.  In  order 
to  confirm  the  conjecture  that  the  high-frequency/low-frequency  (HF/LF)  discriminant 
is  a  depth  of  burial  effect  will  require  calculations  for  greater  overburial  and  smaller 
yields.  Those  explosions  that  Taylor,  et  al.  (1988)  found  to  fail  the  discriminant, 
were  low-yield  explosions  overburied  by  100%  or  more.  Calculations  presented  in  this 
report  indicate  that  both  a  low-frequency  enhancement  is  seen  in  the  normal  buried 
shot  relative  to  the  overburied  shot,  but  that  high  frequency  enhancement  is  also  evi¬ 
dent  in  the  overburied  source.  Calculations  should  be  repeated  to  simulate  smaller 
explosions  with  substantial  overburial  in  a  refined  tuff  model. 
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The  extent  of  the  nonlinear  attenuation  that  affects  the  radiation  of  "pP"  and  SV 
from  the  nonlinear  simulations  should  be  addressed.  Clearly,  the  tuff  models  can  be 
improved  by  introduction  of  additional  layering  to  simulate  the  low  yield  strength 
interbedded  tuff  and  ash  fall  layers  in  the  geologic  sections.  Analysis  indicates  that 
spall  should  be  an  important  contributor  to  regional  seismograms  and  in  particular 
could  be  an  effective  Lg  radiator.  Introduction  of  lower  yield  strength  tuffs  above  the 
working  point  should  increase  the  predicted  explosion  moments  and  increase  the  depth 
to  which  spall  occurs  in  the  simulations.  This  will  increase  the  spall  momentum  and 
increase  the  CLVD  component  of  the  seismic  source.  Additionally,  the  introduction  of 
thin  weak  layers  into  the  model  may  encourage  multiple  spall  separations  and  increase 
SV  radiation. 

Competing  linear  scattering  mechanisms  for  the  generation  of  Lg  by  explosions 
should  be  investigated.  It  is  evident  that  the  biggest  differences  in  the  2D  and  ID  syn¬ 
thetics  are  present  in  the  lowest  group  velocity  modes  of  propagation.  Since  the  low 
velocity  near  surface  layers  of  the  crust  are  responsible  for  these  modes,  it  is  unclear 
where  this  energy  ultimately  goes  since  these  surface  layers  are  rarely  continuous  and 
highly  heterogeneous.  Therefore,  scattering  must  play  an  important  part  in  the  energy 
budget  of  these  modes  and  the  ultimate  fate  of  this  energy  radiated  by  the  shallow 
explosion.  It  has  been  proposed  that  this  energy  is  coupled  into  the  Lg  wavetrain. 

Care  should  be  taken  in  assuming  that  geometrical  spreading  for  Pn  (and  probably 
Sn)  is  frequency  independent  at  high  frequencies.  The  velocity  gradient  (in  any  region 
of  the  world)  below  the  Moho  is  poorly  constrained  by  either  refraction  or  reflection 
data  and  therefore,  only  poorly  known.  The  sensitivity  of  the  frequency  dependence  of 
the  Pn  apparent  attenuation  to  this  gradient  implies  that  strong  trade-offs  can  exist 
between  source  models  and  attenuation  models  for  Pn  unless  the  source  models  are 
independently  constrained.  The  spectral  shape  is  insufficient  to  constrain  the  frequency 
dependence  of  Q  and  spatial  decay  information  is  required  at  multiple  frequencies. 

The  Shagan  River  simulations  discussed  in  this  report  demonstrate  that  explosions 
in  a  simple  high  velocity  brittle  material  may  exhibit  substantial  differences  at  regional 
distances  from  a  point  source  model.  However,  it  is  clear  that  this  simple  model  for 
the  Eastern  Kazakhstan  test  site  can  and  should  be  improved  in  the  light  of  new  infor¬ 
mation  about  the  testing  locality.  Information  about  rock  types,  velocities,  layering 
information,  and  testing  depths  should  be  incorporated  into  the  nonlinear  explosion 
models  for  comparison  with  ID  models  for  explosions  and  earthquakes. 
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Since,  most  magnitudes  or  other  yield  estimators  for  the  Shagan  River  test  site  are 
based  on  m^(Pmax),  m^(Lg),  or  Mj,  we  have  to  relate  the  mt(Pa)  and  ^^(Pb) 
estimates  of  Section  4.2.2  to  other  measures  that  have  been  made  for  a  great  deal  more 
events.  It  is  generally  thought  that  if  various  independent  seismic  "magnitudes"  are 
used,  then  better  yield  estimates  result  due  to  the  cancellation  of  random  errors  on  an 
event  by  event  basis.  We  do  this  by  forming  a  general  linear  model  for  the  estimation 
of  a  "GLM"  unified  magnitude.  We  assume  that  the  slope  between  Mg  and  Mq  or  Mj 
with  body-wave  magnitude  is  0.9,  and  we  use  the  relationship,  -  ^og(^o)  ‘  11-86, 
from  Stevens  (1986)  to  make  a  first  approximation  of  converting  moments  to  Ms-  Mq 
are  taken  firom  Stevens  (1986),  Mj  are  taken  from  Stevens  (personal  communication) 
and  Given  and  Mellman  (1986).  Ms  are  taken  from  from  Sykes  and  Ruggi  (1986). 
We  also  use  the  mi,(Lg)  estimates  of  Nuttli  (1986),  the  P  coda  estimates  of  Gupta 
(1986),  and  the  P  coda  and  rms  Lg  estimates  of  Ringdal  (1987).  Gupta’s  P  coda 
measurements  were  corrected  for  an  estimated  source  spectrum  and  are  theoretically 
proportional  to  moment,  therefore  we  correct  them  to  assuming  a  0.9  slope.  We 
use  two  sets  of  estimates.  12  events  (less  the  cratering  event)  for  which  McLaugh¬ 
lin,  et  al.  (1986)  estimated  maximum  likelihood  WWSSN  magnitudes  for  m^(a), 
mft(b),  and  /Mi,(max),  and  the  ISC  estimates  of  Marshall,  et  al.  (1984).  Since 
fewer  than  15%  of  the  ISC  stations  used  by  Marshall  et  al.  are  WWSSN  stations,  the 
two  sources  of  are  largely  based  on  independent  data  sets.  We  expect  that  the 
nti,  (Pa),  nth  (Pb)  and  (Pmax)  estimates  are  all  highly  correlated. 

In  order  to  use  the  general  linear  model,  we  ignore  correlation  that  surely  exists 
between  the  various  P  measures,  Lg  measures,  and  LR  measures  of  yield.  Since  the 
correlations  are  unknown  at  this  time,  and  the  full  general  linear  model  with  correla¬ 
tion  is  numerically  intractable,  we  are  forced  at  this  stage  to  assume  independence. 
Some  care  has  been  taken  to  balance  the  data  sets  so  that  this  should  not  be  an  overly 
important  consideration.  Given  this  consideration,  we  are  probably  underestimating  the 
standard  errors  of  the  correction  terms  and  the  event  GLM  magnitude. 

The  results  are  shown  below  in  Tables  A.l  and  A.2.  The  baseline  corrections  for 
the  different  source  measures  have  been  adjusted  so  that  the  m^flSC)  estimates  of 
Marshall,  et  al.  (1984)  have  no  net  correction.  For  comparison  with  another  set  of 
magnitudes,  the  average  difference  between  these  magnitude  estimates  and  those  of 
Barker  and  Murphy  (1988)  is  0.00(0.03)  with  an  rms  value  of  0.11. 

From  this  data,  we  can  now  estimate  the  GLM  magnitude  for  the  January  15, 
1965  event  based  on  the  available  mj(Pa),  m^(Pb),  and  mfc(Lg)  estimates.  We 
have  estimated  uncertainties  for  the  m|,(Lg)  of  0.05,  and  Ms  of  0.07.  Stevens  and 
McLaughlin  (1988)  have  estimated  moments  for  the  NTS  cratering  shots  SEDAN  and 
SCH(X)NER.  They  found  that  these  events  had  log(Afo)’s  0.6  log  units  lower  than 
comparable  contained  shots  of  the  same  yield.  We  estimate  an  additional  uncertainty 
in  this  correction  of  0.1  log  units  because  of  various  factors  concerning  the  differences 
between  explosion  coupling  in  different  media.  Therefore,  the  Sykes  and  Ruggi  Ms 
for  the  January  15,  1965  event  should  be  corrected  to  3.9(0.07)  +  0.6(0. 1).  Nuttli’ s 
data  for  SEDAN  show  no  evidence  that  cratering  made  a  difference  in  mi,(Lg).  The 
results  are  summarized  as  follows. 
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nth  (Pa  )=  5.48(0.05)-0. 128-K).661 (0.034)=6.013(0.060) 
nth  (Pb  )=  5.72(0.05)-K).000+0.312(0.030)=6.032(0.058) 
Ms=  0.9*  (3.9(0.07>K).06(0.1)H2.262(0.047)=6.31(0.12) 
nth  (Lg  )=  5.87(0.05>f0.089(0.034)=5.959(0.060) 


We  are  faced  with  forming  some  sort  of  weighted  average  of  the  estimates,  based  on 
their  reliability.  If  we  include  the  Ms  estimate  in  a  weighted  average  then  we  have 

mi, (125kt)  =  6.025(0.03).  The  variance  of  the  Ms  estimate  is  probably  under 
estimated  as  it  is  the  most  uncertain  data  point.  Given  that  we  prefer  the  weighted 
average  of  the  P  and  Lg  estimates,  m^(125kt)  =  6.00(0.03).  Assuming  an 

mj,  :Log(yield)  slope  of  0.9,  this  infers  a  GLM  unified  mfc(150kt)  =  6.07(0.03). 

It  should  be  noted  that  systematic  variations  in  magnitude  differences  (see  Barker 
and  Murphy  1988;  Ringdal,  1987)  have  been  observed  across  the  Shagan  River  test 
site,  and  that  the  January  15,  1965  event  was  located  in  the  southeast  comer  of  the  test 
site  (Rodean,  1979;  Marshall,  et  al.  1984).  This  introduces  an  additional  systematic 
uncertainty  to  the  estimated  150  kt  level  due  to  spatial  variability.  These  rms  varia¬ 
tions  are  typically  about  0.1  magnitude  units  across  the  test  site  based  on  the  rms  error 
in  the  GLM.  Therefore,  6.07(0.11)  is  a  conservative  estimate  of  the  150  kt  GLM 
unified  magnitude  and  it’s  uncertainty.  As  can  be  seen  in  Table  A.2,  a  number  of  the 
largest  events  are  within  the  97.7%  upper  bound  (two  sigma)  of  the  inferred  150  kt 
limit. 

We  have  used  the  nonlinear  calculations  to  infer  ntf,  (a)  and  nth  (b)  corrections  for 
the  January  15,  1965  presumed  cratering  explosion  at  the  Shagan  River  test  site.  In 
order  to  relate  these  corrected  magnitudes  to  more  recent  events  of  interest,  we  have 
developed  a  GLM  unified  magnitude  based  on  m^,,  mfc(Lg),  P  coda,  Ms,  Mq,  and  Mj. 
The  results  are  subject  to  two  sources  of  systematic  error.  The  focusing  defocusing  of 
nth  signals  at  Shagan  could  be  responsible  for  a  baseline  error  in  the  cratering  event 
mj,.  We  have  estimated  this  uncertainity  to  lie  in  the  0.11  magnitude  range.  Finally, 
the  corrections  to  the  m{,(a)  and  m^(b)  magnitudes  derived  from  the  nonlinear  calcula¬ 
tions  are  subject  to  model  assumptions  and  our  knowledge  of  the  geophysical  parame¬ 
ters  of  the  Shagan  River  test  site.  The  fact  that  the  corrections  imply  the  correct 
differences  between  nth  (a)  and  ntf,  (b)  for  contained  and  cratering  events  of  comparable 
magnitude  suggest  that  the  simulations  are  properly  modeling  the  early  portion  of  the 
waveform. 
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TABLE  A.l  GLM  MAGNITUDE  CORRECTIONS 


MAGNITUDE  INTERCEPT  a 

SLOPE 

#EVENTS 

mt(ISC) 

0.0 

0.024 

1.0 

13 

nij,  (Pmax) 

+0.084 

0.018 

1.0 

11 

mt(Pb) 

+0.312 

0.018 

1.0 

11 

nth  (Pa) 

+0.661 

0.022 

1.0 

11 

nib  (Lg) 

+0.089 

0.024 

1.0 

15 

rrihihgl) 

+0.103 

0.022 

1.0 

14 

P  coda 

+4.589 

0.024 

0.9 

7 

P  codal 

+0.003 

0.018 

1.0 

15 

Ms 

+2.262 

0.040 

0.9 

8 

Mq 

+2.476 

0.034 

0.9 

5 

Ml 

+2.238 

0.037 

0.9 

13 

TABLE  A.2  SHAGAN 

RIVER  GLM-MAGNITUDES 

EVENT 

YY/MM/DD 

GLM-MAG 

a 

#MAGS 

1 

68/07/19 

5.355 

0.038 

5 

2 

69/11/30 

6.017 

0.031 

5 

3 

73/07/23 

6.303 

0.023 

9 

4 

73/12/14 

5.908 

0.030 

8 

5 

75/04/27 

5.625 

0.019 

8 

6 

76/07/04 

5.932 

0.021 

10 

7 

76/12/07 

5.777 

0.039 

6 

8 

78/06/11 

5.963 

0.024 

11 

9 

78/09/15 

5.960 

0.020 

11 

10 

79/06/23 

6.160 

0.026 

11 

11 

80/09/14 

6.178 

0.025 

10 

12 

81/09/13 

6.179 

0.033 

5 

13 

81/12/27 

6.027 

0.081 

5 

14 

83/06/12 

6.182 

0.043 

3 

15 

83/10/26 

6.161 

0.037 

4 

16 

84/07/14 

6.193 

0.044 

4 

17 

84/10/27 

6.201 

0.033 

4 

18 

84/12/16 

6.176 

0.037 

4 
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TABT  H  B.la  Frequency  Dependent  Q  Model  SS 


Phase 

DOB 

TOA=l(F 

TOA=20^ 

TOA=3(f 

Log  {Pa) 

200 

0.402 

0.389 

0.370 

LogiPb) 

200 

0.570 

0.553 

0.531 

Log  {Pc) 

200 

0.591 

0.553 

0.519 

200 

0.168 

0.165 

0.161 

200 

0.189 

0.164 

0.149 

Log  {Pa) 

680 

0.600 

0.595 

0.580 

Log{Pb) 

680 

0.791 

0.788 

0.779 

Log  {Pc) 

680 

0.814 

0.802 

0.794 

680 

0.191 

0.192 

0.199 

Pc 

680 

0.214 

0.206 

0.214 
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TABLE  B.lb  Frequency  Dependent  Q  Model  TS 


Phase 

DOB 

TOA=10" 

TOA=20" 

TOA=30® 

Log  {Pa) 

200 

0.043 

0.029 

0.015 

Log{Pb) 

200 

0.178 

0.165 

0.150 

Log  {Pc) 

200 

0.204 

0.176 

0.154 

Log(—) 

200 

0.135 

0.136 

0.135 

Pc 

200 

0.161 

0.148 

0.139 

Log  {Pa) 

680 

0.236 

0.233 

0.221 

Log{Pb) 

680 

0.414 

0.414 

0.408 

Log  {Pc) 

680 

0.455 

0.448 

0.443 

j 

680 

0.178 

0.181 

0.188 

Pc 


680 


0.219 


0.215 


0.223 
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TABLE  B.lc  Frequency  Dependent  Q  Model  TT 


Phase 

DOB 

TOA=10^ 

TOA=20" 

TOA=30" 

Log  (Pa) 

200 

-0.266 

-0.281 

-0.295 

LogiPb) 

200 

-0.151 

-0.160 

-0.169 

Log  (Pc) 

200 

-0.134 

-0.159 

-0.183 

Lcsiy^) 

200 

0.115 

0.121 

0.126 

Pc 

200 

0.132 

0.122 

0.112 

Log  (Pa) 

680 

-0.080 

-0.083 

-0.089 

Log(Pb) 

680 

0.101 

0.100 

0.096 

Log  (Pc) 

680 

0.117 

0.118 

0.132 

T  tPb\ 

Log(-) 

680 

0.181 

0.182 

0.185 

I  \ 

680 

0.197 

0.200 

0.221 
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TABLE  B.ld  Frequency  Independent  Q  Model  ( t*  =  0.5  ) 


Phase 

DOB 

TOA=10® 

TOA=20® 

TOA=30® 

LogiPa) 

200 

0.047 

0.031 

0.017 

LogiPb) 

200 

0.196 

0.185 

0.174 

Log  {Pc) 

200 

0.239 

0.216 

0.198 

t 

U>gij^) 

200 

0.151 

0.154 

0.156 

Log(-^) 

200 

0.194 

0.185 

0.180 

Log  {Pa) 

680 

0.234 

0.231 

0.226 

Log{Pb) 

680 

0.440 

0.440 

0.437 

Log  {Pc) 

680 

0.498 

0.494 

0.497 

Log{-) 

680 

0.206 

0.208 

0.211 

r  /P^  \ 

680 

0.214 

0.263 

0.271 
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TABLE  B.le  Frequency  Independent  Q  Model  ( f*  =  0.7  ) 


Phase 

DOB 

TOA=10" 

TOA=2(r 

TOA=30" 

Log  {Pa) 

200 

-0.237 

-0.251 

-0.265 

LogiPb) 

200 

-.097 

-0.105 

-0.115 

Log  (Pc) 

200 

-0.065 

-0.091 

-0.119 

r  Pbs 

200 

0.139 

0.146 

0.150 

Pc 

200 

0.172 

0.160 

0.146 

Log  {Pa ) 

680 

-0.055 

-0.058 

-0.064 

Log{Pb) 

680 

0.153 

0.151 

0.147 

Log  {Pc) 

680 

0.182 

0.190 

0.210 

I 

Logi.J^) 

680 

0.208 

0.209 

0.210 

Pc 

680 

0.237 

0.247 

0.273 

Unfllt«r*d  2d  va  Id  aolutlona  at  toaa39.  S  dab 
Data:  11/  2/19M  Tlmmt  l«tSa:63 


Radiated  displacement  P  waveforms  for  the  three  2D  Pahute  Mesa  tuff  simulati 
at  30  takeoff  angle  compared  to  the  ID  elastic  source  model. 


Figure  B.3.  Same  as  B.2  except,  autoscaled  to  “b”  phase. 
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APPLICATION  OF  THE  REPRESENTATION  THEOREM  AND  USE  OF  THE 
“RING-LOAD”  GREEN’S  FUNCTION  FOR  REGIONAL  SYNTHETIC  CALCULATIONS 
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By  use  of  the  representation  theorem,  the  results  of  finite  difference  simulations 
may  be  coupled  to  Green’s  functions  for  regional  and  teleseismic  distances  (Figure 
C.l).  The  method  has  been  used  by  Day,  et  al.  (1986)  with  modal  suirunation 
Green’s  functions  to  study  surface  waves  from  nonlinear  finite  difference  simulations, 
and  recently  by  Day,  et  al.  (1986b)  to  calculate  teleseismic  P  waves  from  finite 
difference  source  simulations  using  teleseismic  layered  stack  propagators.  Regional 
waves  are  too  complicated  to  be  modeled  properly  with  modal  summation  procedures 
or  with  simple  far-field  Haskell-Thompson  propagators. 

In  practice,  we  use  Green’s  functions  for  a  layered  Earth  computed  with  a 
wavenumber  integration  algorithm  (see  Luco  and  Apsel,  1983),  but  the  method  applies 
in  general  to  (jreen’s  functions  computed  by  other  methods.  In  this  way,  we  simulate 
complicated  nonlinear  as  well  as  linear  wave  propagation  in  a  grid  using  finite 
differences  and  then  propagate  the  wavefields  to  large  distances  in  a  layered  Earth. 
The  procedure  has  a  number  of  applications  for  seismic  scattering  problems  as  well  as 
for  use  with  nonlinear  seismic  sources. 

The  representation  theorem  (Aki  and  Richards,  1980)  states  that  the  displacement 
field  at  location  io  roay  be  given  by  an  integral  over  the  displacements,  u(x.),  and  trac¬ 
tions,  Og(i),  over  a  closed  surface,  X.  In  the  frequency  domain,  we  can  write  that, 

“  M{l)  C:Vg,(jj£o)]  (1) 

where  g  (i  zo)  is  the  displacement  at  i  due  to  a  point  force  in  the  Xi  direction  applied 

at  xq,  C  is  the  elastic  tensor,  m(i)  is  the  displacement  at  £  from  the  finite  difference 
simulation,  and  is  the  stress  tensor  on  X  from  the  finite  difference  simulation.  In 
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Figure  C.l.  2D  nonlinear  finite  difference  calculations  are  used  to  compute  displacements  and 

tractions  on  a  monitoring  surface.  The  motions  and  forces  on  the  monitoring 
surface  are  then  used  to  propagate  both  regional  and  teleseismic  far-field  body 
waves  from  the  2D  nonlinear  sources  using  the  elastodynamic  representation 
theorem. 
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cylindrical  coordinates  the  integrand  becomes, 

r(i io)  =  gir'^r -  M :7gi  (2) 

and  equation  (1)  becomes 

Ui(js.o)  ==  fdS^r(xJS.o)  =  j^rdrdz d^^(iJ^o)  (3) 

where  =  f'C  n,  Xj  =  i  c  n,  and 

M  =  nU‘C  =  X/n-R  +  ^.(ni<  +  un). 

The  integration  is  performed  over  a  cylindrical  surface  composed  of  a  surface  of  con¬ 
stant  radius,  r,  referred  to  as  the  "side",  and  a  surface  of  constant  depth,  z,  referred  to 
as  the  "bottom".  For  the  side, 

M  =  -  Ur (Xf  +  liiff)  -  +  if), 

and  for  the  bottom. 


Af  =  -  M,(X/  +  2^f£)  -  [uirifi  +  if). 

The  gradient  of  the  Green’s  function  in  cylindrical  coordinates  becomes, 

V*,  = +  W(-^)  +  «'’(%-)  +  «(%■)  +  **(%•)•  f"*) 

By  reciprocity,  gir(XJLo)  can  be  re-interpreted  as  the  i’th  component  of  displace¬ 


ment  at  2o  <lue  to  the  point  source  in  the  £  direction  at  2.  To  compute  the  azimuthal 


integral  in  equation  (3),  note  that  the  displacements  and  tractions  in  the  integrand  are 
independent  of  0  (axi-symmetric  finite  difference  output)  and  the  integral  may  be  po*- 
formed  over  the  Green’s  functions  and  the  derivatives  of  the  Green’s  functions  with 


respect  to  depth  and  range.  These  functions  depend  on  ^  by  virtue  of  their  dependence 
on  2o  (tfic  receiver  location)  and  we  will  label  the  azimuthally  integrated  Green’s  func¬ 
tion  terms  as  Gjo  and  dpG,o  (a  =  r,z,  and  P  =  r,z).  Gi„(r,,z,,ro,zo)  is  interpreted  as 
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the  i’th  displacement  at  (rQ,zo)  due  to  the  ring  loads  /a(r,z)  =  8(r-r,)8(z-z,)  in  the 
a-f  and  z  directions. 

Putting  equations  (2),  (3),  and  (4)  together  using  the  azimuthally  integrated 
Green’s  functions,  then  the  contribution  from  the  side  at  radius  r=a  becomes, 

Ui  (r  o,z  o)  =  «  fi^iz  "Cz  Giz  '^z  (5a) 

+  aj[(X  +  2|i)u^9,G,y  +  +  djGiy)]dz 

and  from  the  bottom  at  depth  z, 

Ui  {rQ,z  o)  =  J[G,>  +  G^  T,  ]rdr  (5b) 

+  jlXu^d^Gi,  +  |XM,(8^G,y  +  9^0;^)  +  u^((X  +  2pL)dgGig  +  ■yGiy)]rdr 

Monitored  displacements  and  tractions  on  the  cylindrical  surface  are  saved  at 
individual  time  steps  of  the  finite  difference  calculation  for  a  plane  of  constant  depth 
and  a  set  of  grid  locations  of  constant  radius.  The  Green’s  functions  and  their  deriva¬ 
tives  are  computed  for  ring  loads  using  a  modified  version  of  a  wavenumber  integra¬ 
tion  algorithm.  The  far-field  synthetics  are  then  formed  from  the  sum  of  the  appropri¬ 
ate  convolutions.  The  wave-number  integration  program  has  been  ported  to  a  vector 
processor  computer  and  substantial  improvements  have  been  made  since  many  of  the 
source  locations  at  constant  depth  share  the  same  wavenumber  layer  matrix  and  only 
differ  in  range. 

The  results  have  been  tested  against  a  half-space  analytic  solution  of  a  point 
explosion  source  observed  at  various  take-off  angles  in  the  half-space  as  well  as  by  use 
of  the  far-field  continuation  procedure  used  in  Day,  et  al.  (1986)  and  used  in  Section 
IV  of  this  report.  Tests  have  also  been  done  by  comparison  with  a  similar  program 
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that  computes  the  excitation  of  Rayleigh  modes  and  propagates  these  to  the  far-field 
using  the  output  of  the  finite  difference  calculations.  Additional  tests  include  the  sys¬ 
tematic  removal  of  monitored  finite  difference  grid  locations  in  the  summation  pro¬ 
cedure.  As  should  be  expected,  decimation  produces  the  effect  that  the  proper  summa¬ 
tion  begins  to  break  down  at  high  frequencies  such  that  the  wavefield  becomes  spa¬ 
tially  aliased. 
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